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SUMMARY 


This  report  presents  a  general  method  for  predicting  and  bounding 
casualties  from  both  prompt  and  delayed  effects  produced  by  nuclear  attacks 
against  population.  We  summarize  the  major  findings  and  recommendations 
as  follows: 

1.  The  feasibility  of  interrelating  analytically  the  major  ingredients 
of  damage/ casualty  prediction  has  been  demonstrated  by  deriving  tractable, 
explicit  analytic  representations  of  the  interplay  among  the  characteristics 
of  the  attack,  targets,  environment  and  outcome. 

2.  Our  approach  for  determining  the  effective  blast  protection  of  a 
population  with  mixed,  non-uniform  protection  produces  stable,  consistent 
values  of  protection  hardness  and  accurate  estimates  of  casualties. 

3.  Relocation  schemes  which  evacuate  a  region's  population  to  locations 
where  at  least  some  people  already  reside  have  been  characterized  abstractly 
and  analyzed  parametrically;  in  particular,  the  conditions  for,  and  con¬ 
sequences  of,  self-contained  regional  relocation  have  been  determined. 

4.  It  is  highly  unlikely  that  any  realistic  approach  could  relocate 

the  bulk  of  the  nation's  population  at  population  density  less  than  1,000  to 
2 

1,500  people/mi  ,  i.e.,  some  2  to  3  times  the  population  density  corresponding 
to  spreading  the  entire  population  uniformly  over  the  entire  populated  area. 

5.  To  achieve  moderate  survivor  fractions  for  attacks  against  popula¬ 
tion  with  several  thousand  blast-equivalent  megatons,  significant  blast 


hardening  is  required,  even  with  the  most  defense-optimistic  assumptions 
about  relocation. 

6.  For  attacks  against  population,  the  influence  of  delayed-effects 
hardness  is  generally  weak,  compared  to  the  influence  of  prompt-effects 
hardness  and  dispersion. 

7.  In  view  of  the  results  and  conclusions  presented  above,  we  recommend 

that 

a)  the  demographics  of  residential  population  be  analyzed  to 
explore  the  validity  of  the  apparent  universality  of  normalized  rank-ordered 
population  distribution  for  different  geographic  regions  and  to  provide  a 
basis  for  interpre  .ig  the  rank-ordered  population  distribution  in  other 
commonly  used  formats  (e.g.,  in  terms  of  political  subdivisions) 

b)  the  general  analytic  approach  and  results  be  extended,  refined, 
and  applied,  e.g.,  by  considering  attacks  more  general  than  attacks  against 
population,  tighter  bounds  on  combined-effects  casualties,  and  specific  candi¬ 
date  relocation  schemes. 


1.  INTRODUCTION 

This  report  completes  the  presentation  of  results  obtained  by  La  Jolla 
Institute  in  a  continuing  investigation  supported  by  Contracts  No.  DCPA01- 
73-C-0183  and  DCPA01 -79-C-0244.  Earlier  results  are  presented  elsewhere  [Ref.  1] 
and  are  summarized  below  as  needed.  The  overriding  objective  of  this  effort 
has  been  to  interrelate  analytically  the  major  ingredients  of  damage/casualty 
prediction  by  deriving  tractable,  explicit  analytic  expressions  of  the  inter¬ 
play  among  the  characteristics  of  the  attack,  targets,  environment  and  outcome. 

In  working  toward  this  objective,  we  have  put  special  emphasis  upon  nuclear 
attacks  against  population  and  upon  the  civil  defense  measures  of  relocation 
and  hardening.  The  spirit  of  the  investigation  has  been  to  seek  not  specific 
numerical  results  but  rather  tractable  analytic  expressions  and  thereby  to 
understand  parametric  trends  and  to  explore  the  meaningful  ness ,  accuracy,  and 
sensitivity  of  the  results.  We  have  accordingly  sought  solutions  for  not  only 
the  direct  problem  (i.e.,  predicting  the  outcome  as  a  function  of  the  parameters 
characterizing  the  target,  attack,  and  environment)  but  also  for  the  inverse 
problem  of  expressing  any  ingredient  of  the  problem  as  a  function  of  the  others. 
In  particular,  we  have  sought  results  which  prescribe  the  hardening/dispersion 
combinations  sufficient  to  guarantee  a  specified  outcome.  Such  results  clearly 
provide  a  rational  basis  for  designing  passive  defense  measures  and  for  balancing 
hardness  and  dispersion  measures. 
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In  Section  2  below,  we  begin  by  summarizing  our  general  analytic  results 
(including  solutions  of  the  direct  and  inverse  problems  discussed  above)  and 
their  conceptual  foundations.  We  then  discuss  a  major  element  of  each  of  the 
two  generic  civil  defense  techniques  of  hardening  and  dispersion;  speeifical ly, 
we  present  results  for  characterizing  effective  blast  hardness  and  relocation 
schemes  which  evacuate  people  to  locations  where  at  least  some  people  already 
reside.  In  closing  Section  2,  we  illustrate  the  application  of  our  general 
analytic  results  by  presenting  a  number  of  specific  examples  for  prompt  and 
delayed  casualties  and  for  in-place  and  relocated  populations.  In  Section  3, 
we  draw  conclusions  and  make  recommendations. 
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2.  STRUCTURE  OF  THE  DAMAGE/ CASUALTY  PREDICTION  PROBLEM 

In  this  section,  we  summarize  and  illustrate  our  results  concerning  the 
structure  of  the  damage/casualty  prediction  problem. 

2.1  Summary  of  Analytic  Results 

Unless  otherwise  noted,  the  results  presented  in  this  subsection 
are  derived  in  Ref.  1,  frequently  in  considerably  greater  generality  than 
indicated  below;  for  complete  derivations  of  the  results  (including  the  range 
of  their  validity),  this  reference  should  be  consulted.  We  begin  by  recapitu¬ 
lating  the  ingredients  of  the  problem. 

For  an  attack  against  population,  the  target  is  characterized  by 
its  joint  population/hardness  distribution.  The  hardness  is  specified  by 
four  parameters:  the  prompt-effects  hardness  by  a  single  effective  damage 
overpressure,  Ap$  5Q;  the  delayed-effects  (i.e.,  fallout)  hardness  by  the 
effective  protection  factor  <p>  (i.e.,  the  attenuation  of  the  free-field  dose 
by  an  individual's  surroundings)  and  by  the  mean,  m,  and  standard  devia¬ 
tion,  ap,  defining  the  "universal"  Gaussian  distribution  of  damage  probability 
vs.  total  equivalent  biological  dose  to  which  an  individual  is  exposed.  (It 
will  be  recalled  that  the  pair  [Dgff  m,aD]  is  customarily  taken  to  be 
[450R,  75R]  and  [250R,  50R]  for  fatalities  and  casualties,  respectively.)  The 
spatial  distribution  of  population  is  specified  by  the  distribution,  F,  of 
cumulative  population,  si,  versus  cumulative  area.  A,  obtained  by  successively 
subdividing  the  target  region  into  rectangular  cells,  noting  the  population 
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and  area  of  each  cell,  rank-ordering  the  cells  by  population  density  (i.e., 
by  population  divided  by  area),  and  summing  over  the  rank-ordered  cells.  Thus 
the  population  distribution  (normalized  by  the  total  regional  population,  n  t) 
is  given  by 


0) 


where  A*  is  a  suitable  normalizing  area.  The  population  density,  u,  is  given  by 


03  = 


dfl  ntot  P, /A, 
dA  "  A*  *■  lA*; 


(2) 


From  the  construction  of  F,  it  is  clear  that  the  population  density  (and  hence 
F* )  is  positive  and  decreasing,  so  that  the  distribution  F  is  convex.  In  sum, 
then,  the  target  is  characterized  by  the  four  parameters  specifing  the  hardness, 
plus  the  two  parameters  and  A*  which,  along  with  any  parameters  character¬ 
izing  the  dimensionless  distribution  F,  specify  the  spatial  distribution  of  the 
population. 

An  attack  against  population  targets  the  first  weapon  against  the  most 
populous  cell  and  sequentially  targets  cells  rank-ordered  by  population  density. 

To  characterize  an  attack  against  population,  it  is  clearly  necessary  to 
specify  the  total  number,  N,  and  the  burst  height  of  the  attacking  weapons 
and  to  specify  the  yield,  Y,  and  fission  fraction,  F/Y,  of  the  individual  weapons. 
(To  avoid  confusion  with  the  symbol  introduced  in  Equation  1  to  describe  the 
normalized  population  distribution,  we  refer  to  the  fission-yield  exclusively 
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in  the  fission-fraction  combination  F/Y.)  Since  the  effective  damage  over¬ 
pressure  is  typically  less  than  10  psi  and  never  exceeds  tens  of  psi ,  it  suffices 
to  consider  two  burst-heights:  an  airburst  optimized  for  the  damage  overpressure 
maximizes  the  prompt-effects  damage  but  produces  no  fallout  and  hence  no  delayed- 
effects  damage;  a  surface  burst  maximizes  the  delayed-effects  damage  and  produces 
the  same  prompt-effects  damage  as  any  fallout-producing  detonation.  Therefore, 
to  complete  the  attack  characterization,  it  suffices  to  specify  either  the  number, 
N^,  of  surface  bursts  or  the  number,  N^,  of  optimized  airbursts,  since 

N  «  Na+Ns  (3) 

In  sum,  then,  the  attack  is  characterized  by  four  parameters:  Y,  F/Y,  N,  N<.. 

The  environment  for  prompt-effects  damage  is  characterized  by  the 
ambient  pressure,  PQ,  while  the  environment  for  delayed-effects  damage  is 
characterized  by  the  "effective  fallout"  wind,  W,  and  cross-wind  shear,  Sc, 
where,  here  and  throughout  this  report,  a  symbol  with  underscore  denotes 
a  vector  quantity,  while  the  same  symbol  without  the  underscore  denotes  the 
magnitude  of  the  vector.  In  our  approach,  by  suitably  formulating  faTlout 
phenomenology  and  employing  bounding  arguments,  we  are  able  to  finesse  the 
effects  of  W  and  Sc,  and  to  obtain  delayed-effects  damage  results  independent 
of  these  parameters. 

The  outcome  of  the  attack  is  characterized  by  the  number  of  fatalities, 
or  casualties,  nca$.  In  stating  results,  additional  affixes  have  the 
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obvious  interpretations:  e.g.,  Qfatjprompt  denotes  the  fatalities  from  prompt 
effects  while  com^  denotes  the  fatalities  from  combined  (i.e.,  prompt 
plus  delayed)  effects.  The  outcome  may  equivalently  be  expressed  in  terms 
of  the  survivor  fraction  e  given  by 


« =  1  -  (4) 

“tot 

The  parameters  introduced  above  characterize  the  ingredients  of  the 
damage/casualty  prediction  problem.  For  attacks  against  population  by  con¬ 
temporary  nuclear  weapon  delivery  systems,  we  find  that  these  ingredients  are 
interrelated  via 


fat,promi 


vs  S50,S+(1'vs)s50 


>.*])- F(: 


50, A 
A* 


fat, comb 


[l+4>] 


where 


[2  j  offense  /  i  lower/ 

t  *  j  j  for  j  -  conservative,  i.e.,  j  >  bounds  on  fatalities 

' tt  ’  defense)  'upper' 
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■  "A  W/P,)'"^  ) 


The  parameters  s5Q  A  and  SgQ  $  express  the  distance  to  which  the 

damage  overpressure  Ap$  extends  for  an  optimized  airburst  and  for  a  surface 

burst,  respectively.  The  parameters  y  ,  and  A  arise  from  our  formulation 

max  o 

of  fallout  phenomenology  and,  in  particular,  from  our  result  that,  for  a  single 
detonation,  the  isodose  area  (i.e.,  the  area  within  an  isodose  contour,  or, 
equivalently,  the  area  within  which  a  given  dose  is  exceeded)  corresponding  to 
a  free-field  dose  D  is  given  by 


A .  A  =AJiD\  A  /aD\ 

isodose  \A  /  D  \A  ) 

(17) 

where 

A  =  A0(F/Y) 

(18) 

and  the  normalizing  area  *  is  given  by 

■  A  =  ouS  WT2 

H  c 

(19) 

and  the  parameters  cH  and  T  depend  on  yield  only  and  can 

[Ref.  2]  for  yields  in  the  megaton  range  as 

be  well 

approximated 

0|_|  *  [7.5  +  1.5  log^g  Y]kft 

Y  in  Mt 

(20) 

T  =  [7.5  +  1.66  log1Q  Y]hr 

Y  in  Mt 

(21) 
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and  the  functions  4  and  v  are  shown  in  Figure  1.  As  may  be  seen  from  Figure  1, 
the  parameter  of  Equation  14  is  indeed  the  maximum  value  of  the  function  y,  i.e., 

v  =  °-4869  =  y(w=w*=0.6834)  (22) 

By  selecting  aim-points  at  (or  above)  the  center  of  the  N  most  densely 
populated  cells,  an  attack  with  N  weapons  is  guaranteed  to  achieve  the  prompt- 
effects  fatalities  given  by  Equation  5,  with  t=2.  Conversely,  an  attack  with 
N  weapons  cannot  achieve  greater  prompt-effects  fatalities  than  those  given  by 
Equation  5,  with  t=tt.  Therefore,  the  actual  number  of  prompt-effects  fatalities 
is  given  by  Equation  5  for  some  value  of  t  between  2  and  it;  moreover,  the 
values  of  2  and  ir  for  t  provide  the  bounds  indicated  by  Equation  7.  It  is 
clear  that  an  attack  employing  airbursts  exclusively  (i.e.,  v$=0)  is  feasible 
and  maximizes  the  prompt-effects  fatalities.  The  upper-bound  on  combined  fatali¬ 
ties  as  expressed  by  Equation  6  is  derived  by  upper  bounding  the  size  and  pop¬ 
ulation  of  the  fallout  region;  although  lower  bounds  or  alternative  upper  bounds 
could  also  be  derived,  the  result  expressed  by  Equation  6  is  especially  powerful, 
as  we  describe  below;  here,  and  throughout  this  report,  we  follow  the  convention 
that  "Equation  is  a  name  for  a  mathematical  relationship  which  may,  but  need 
not  be,  an  equality. 

It  is  clear  that  by  presenting  the  outcome  as  an  explicit  analytic 
function  of  the  problem's  other  ingredients,  Equations  5  and  6  solve  the  "direct" 
problem  discussed  above.  Moreover,  the  solution  of  the  inverse  problem  is 
readily  inferred.  Specifically,  Equation  5  immediately  implies  the  equivalent 
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2/3 

blast  yield  NY  '  (popularly  described  as  "equivalent  megatons")  needed  to 
achieve  a  specific  outcome  against  a  specific  target.  Likewise,  the  combina¬ 
tion  of  dispersion  and  prompt-effects  hardness  to  achieve  a  specified  outcome 
is  given  by 

<»> 

where 

Acum(2)  =  A*F-1(z)  »  most  densely  populated  area 

containing  fraction  z  of  total  popu¬ 
lation,  0<z<l  (24) 

Clearly  any  joint  hardness/dispersion  distribution  satisfying  Equation  23 
guarantees  the  specified  outcome  (i.e.,  prompt-effects  survivors)  against  the 
specified  attack.  If  either  the  hardness  or  the  spatial  distribution  of  the 
population  is  specified  independently,  then  the  other  parameter  is  readily 
inferred  from  Equation  23. 

From  Equation  6,  it  is  clear  that  the  combined  (prompt  plus  delayed) 
fatalities  from  a  mixed  surface/airburst  attack  cannot  exceed  the  combined 
(exclusively  prompt)  fatalities  from  a  pure  airburst  attack  if  <p  .is  negative; 
from  Equations  11  and  13,  the  condition  that  41  be  negative  is  equivalent  to 
the  condition 


>  ;  =  _ maxAo/Deff  (25) 

F/Y  ^  /.  2  2. 

T  S50,A  '  s50,S) 

From  Equations  9,  10  and  15  (and  the  fact  that  the  denominator  in  Equation  15 
exhibits  negligible  variation  with  yield),  it  is  clear  that  the  parameter 
depends  only  on  ap$  ^  and  Y  with  the  yield  dependence  given  by 

i*  *  Y1/3  (26) 

3 

Since  the  yield  of  contemporary  strategic  weapons  lies  within  a  factor  of  3 

about  one  megaton  (i.e.,  between  37  kt  and  27  Mt),  varies  with  yield  by 

at  most  a  factor  of  3.  Moreover,  for  attacks  which  produce  the  same  prompt- 

effects  damage  and  which  therefore  have  the  same  blast-equivalent  yield 
2/3 

NY  ,  it  follows  from  Equation  26  that 

ip*  'v  N-^  for  constant  NY^3  (27) 

The  significance  of  Equations  6  and  25  lies  in  the  fact  that  radiation 
shielding  sufficient  to  satisfy  Equation  25  accomplishes  as  much  as  any  amount 
of  shielding  can  accomplish--in  the  sense  that  any  additional  shielding  can 
be  obviated  by  the  offense's  simply  raising  the  burst  height  and  thereby  pro¬ 
ducing  greater  casualties.  The  critical  value  of  shielding  (as  characterized 

1  /3 

by  the  attenuation  <p>  or  protection  factor  PF)  depends  only  on  yield  (as  Y  ) 
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and  on  the  damage  overpressure  and  is  completely  independent  of  all  the  other 
variable  ingredients  of  the  problem--independent  of  how  the  population  is 
distributed  in  an  absolute  sense,  independent  of  how  the  population  is  distributed 
in  a  rank-ordered  sense,  independent  of  the  size  of  the  attack,  independent  of 
where  the  detonations  occur,  independent  of  what  the  winds  and  weather  are, 
independent  of  what  the  number  of  fatalities  is. 

An  approach  to  selecting  performance  criteria  for  passive  defense 
emerges  naturally  from  the  results  presented  above:  first  the  dispersion  and 
prompt-effects  hardness  are  prescribed  by  Equation  23,  then  the  delayed-effects 
hardness  (corresponding  to  the  prescribed  prompt-effects  hardness)  is  given 
by  Equation  25.  In  this  approach,  the  required  dispersion  and  prompt-effects 
hardness  depend  upon  the  desired  outcome  and  upon  the  size  of  the  attack, 
while  the  delayed-effects  hardness  depends  significantly  upon  only  the  prompt- 
effects  hardness.  Delayed-effects  hardness  better  than  that  prescribed  by 
Equation  25  may  reduce  casualties  from  some  attacks,  but  obviously  cannot 
possibly  reduce  casualties  from  the  pure-airburst  attack.  Moreover,  if  the 
delayed-effects  hardness  satisfies  Equation  25,  then  even  if  the  attack  size 
increases  beyond  the  design  level,  it  remains  true  that  no  attack  can  produce 
more  casualties  from  combined  effects  than  those  produced  by  prompt-effects 
from  a  pure  airburst  attack--even  though  the  latter  casualties  increase. 

In  the  results  presented  above,  the  target  hardness  which  can  be  modi¬ 
fied  by  civil  defense  efforts  is  represented  by  two  parameters:  Ap£  ^  and  <o>. 

In  reality,  the  protection  of  the  entire  population  does  not  possess  a  single  homo¬ 
geneous  hardness  for  prompt  or  delayed  effects,  but  can  be  decomposed  into  a 
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finite  number  of  (mutually  exclusive  and  cumulatively  exhaustive)  classes 
of  homogeneous  hardness.  Given  the  latter  more  realistic  and  fine-grained 
distribution  of  hardness,  it  is  possible  (as  shown  in  Ref.  1)  to  define  an 
equivalent  uniform  hardness  so  that  Equations  5  and  6  and  their  implications 
remain  valid.  For  example,  if  there  are  M  exposure  environments  for  delayed 
effects,  with  each  individual  such  environment  characterized  by  the  protec¬ 
tion  factor  p.  and  containing  the  fraction  e.  of  the  total  population,  then 

J  J 

±-t  £ 

j=l  J 

Similarly,  if  each  exposure  environment  for  prompt  effects  contains  fraction 
a.  of  the  total  population  and  is  characterized  by  the  damage  overpressure 

J 

Ap  cn  then  an  analogous  (but  slightly  more  complicated)  expression  gives 
the  equivalent  uniform  hardness  Ap  .n  in  terms  of  all  the  a.  and  Ap  rn  •; 
we  discuss  and  illustrate  this  expression  in  Section  2.2  below. 

Dispersion,  the  other  target  property  which  appears  in  the  results 
presented  above  and  which  can  be  modified  by  civil  defense  efforts,  is  repre¬ 
sented  by  the  rank-ordered  distribution  F  and  its  inverse.  It  is  clear  that 
any  physically  realizable  spatial  distribution  of  population  determines  a 
unique  rank-ordered  F;  the  same  rank-ordered  F  could  (in  principle,  if  not  in 
practice)  arise  from  different  spatial  distributions  of  population.  The 
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distribution  of  residential  (i.e.,  "in-place")  population  is  the  single  most 
important  distribution,  both  for  its  own  sake  and  as  a  starting-point  for 
relocated  distributions.  By  analyzing  the  residential  population  of  the  entire 
nation  and  of  three  sub-regions  (consisting  of  a  12-state  Northeastern  complex, 
a  4-state  Midwestern  complex,  and  the  single  state  of  California,  respectively), 
we  find  [Ref.  1]  that  the  residential  population  of  each  of  these  regions 
determines  a  rank-ordered  distribution  which  is  well  represented  by  F=Fq  where 


F0(w)  ■  V 


-w 


X  =  0.63 


(29) 


The  Northeastern  complex  consists  of  Connecticut,  Delaware,  District  of  Columbia, 
Maine,  Maryland,  Massachusetts,  New  Hampshire,  New  Jersey,  New  York,  Pennsylvania, 
Rhode  Island,  and  Vermont;  the  Midwestern  region  consists  of  Illinois,  Indiana, 
Michigan  and  Ohio.  Moreover,  we  find  that  the  total  populated  area  of  each 
region  is  well  approximated  as  10A*.  Consequently,  different  regions  share  the 
same  normalized  rank-ordered  distribution  of  residential  population  and  differ 
only  in  the  amount  of  unpopulated  area  and,  of  course,  in  the  absolute  scale  of 
the  distribution  as  reflected  in  and  A*.  In  Section  2.3  below,  we  analyze 
the  implications  of  these  results  for  relocation  schemes  which  evacuate  people 
to  locations  where  at  least  some  people  already  reside. 

While  the  preceding  summary  of  analytic  results  is  a  complete  state¬ 
ment  of  the  results  quoted  and  occasionally  suggests  the  derivation  thereof, 
a  more  coherent  presentation  of  their  derivation  may  be  desired.  To  this  end, 
we  now  present  an  overview  (adapted  from  Ref.  1)  of  their  conceptual  founda¬ 
tion.  For  additional  details,  the  reader  should  consult  the  complete  reference 
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which  provides  a  precise  statement  of  the  damage/casualty  prediction  problem 
for  arbitrary  targeting  by  the  attacker  and,  for  attacks  against  population, 
presents  a  general  method  for  predicting  and  bounding  casualties  from  both 
prompt  and  delayed  effects. 

A  precise  and  rigorous  statement  of  the  damage/casualty  prediction 
problem  results  from  adopting  the  conceptual  framework  wherein  targets  are 
susceptible  in  a  finite  number  of  damage  levels  to  a  finite  number  of  different 
"loading"  types  and  are  found  in  a  finite  number  of  mutually  exclusive  exposure 
environments,  with  each  exposure  environment  characterized,  for  each  damage 
level,  by  its  susceptibility  to  each  of  the  loadings.  A  specific  damage/ 
casualty  prediction  technique  is  essentially  defined  by  its  treatment  of  eight 
critical  characteristics.  We  illustrate  this  viewpoint  by  describing  in 
explicit  detail  how  each  of  these  eight  characteristics  is  treated  in  four 
specific  techniques  which  are  in  current  use  and  which  collectively  span  the 
range  of  diversity,  complexity,  and  detail  to  be  found  in  contemporary  approaches. 

To  introduce  the  elements  of  the  original  approach  presented  here 
and  to  lay  the  foundation  for  the  general  method  developed  here,  we  orqanize 
the  discussion  around  the  eight  critical  characteristics  discussed  above. 

These  methodological  auxiliaries  to  the  major  results  show  that 

a)  a  step-function  ("cookie-cutter")  is  a  valid  representation  of 
prompt-effects  damage  probability,  this  result  is  derived  by  theoretical  argu¬ 
ments  and  verified  by  comparison  with  detailed  numerical  results 

b)  damage  degradation  due  to  delivery  system  inaccuracy  is  negli¬ 
gible  for  the  range  of  targets  and  attacks  considered  here;  this  result 
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emerges,  in  the  limit,  from  a  general  analysis  which  demonstrates  that  prob¬ 
ability  of  damage  to  point  and  area  targets  depends  on  at  most  two  over¬ 
pressures  (characterizing  the  offense  and  defense,  respectively)  and  is  tightly 
bounded  by  prescribing  merely  the  ratio  of  these  two  characteristic  over¬ 
pressures 

c)  the  probabilistic  nature  of  casualty  prediction  does  not  intro¬ 
duce  significant  uncertainty  into  the  predictions 

d)  for  fallout-producing  detonations,  the  magnitude  and  location 
of  the  maximum  dose  and  the  variation  thereof  with  wind  can  be  expressed  in 
tractable  analytic  terms,  valid  for  all  physically  realistic  winds;  these 
results  are  derived  from  the  general  WSEG-10  fallout  phenomenology 

e)  the  area  within  isodose  contours  from  fallout-producing  detona¬ 
tions  can  be  expressed  in  simple  analytic  terms  and  upper-bounded  independently 
of  wind  structure;  this  result  follows  from  reformulating  WSEG-10  fallout 
phenomenology  in  a  dimensionless  form  which  provides  "universal"  isodose 
contours 

f)  the  normalized  distribution  of  residential  population  (rank- 
ordered  by  population  density)  versus  cumulative  area  occupied  by  that  popu¬ 
lation  is  essentially  the  same  for  the  regions  analyzed  (i.e.,  the  continental 
United  States,  a  12-state  Northeastern  region,  a  4-four  state  Midwestern  region, 
California)  and  is  described  by  a  simple  analytic  expression;  preliminary 
results  indicate  that  in  all  regions,  the  entire  populated  area  is  approximately 
twenty  times  the  area  containing  the  most  densley  populated  half  of  the  total 
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population  and,  consequently,  that  regions  differ  in  the  amount  of  unpopulated 
area  which  they  contain,  but  do  not  differ  in  the  way  in  which  populated  area 
is  filled. 


Because  different  geographic  regions  appear  to  share  the  same  normal¬ 
ized  rank-ordered  distribution  of  residential  population,  relocation  strategies 
which  evacuate  people  to  locations  where  at  least  some  people  already  reside 
are  more  nearly  similar  than  dissimilar  for  different  regions;  considerations 
of  relocation  on  a  national  scale  may  require  that  different  parametric  ranges 
of  the  common  normalized  distribution  be  emphasized  in  different  geographic 
regi ons . 

The  application  of  our  general  approach  for  predicting  prompt-effects 
casualties  to  40  specific  families  of  attacks  (corresponding  to  every  combina¬ 
tion  of  four  geographic  regions,  two  burst-heights,  and  five  aim-point  spacings) 
against  residential  population  has 

a)  produced  numerical  results  wherein  the  joint  population/hardness 
distribution  is  represented  with  the  most  fine-grained  detail  available  and 
with  the  greatest  defense  optimism  (i.e.,  "best  available  shelter") 

b)  corroborated  analytic  simplifications  in  the  general  approach 
whereby  arbitrary  joint  population/hardness  distributions  are  successively 
reduced  to  the  more  tractable  distributions  of  multiple-uniform-hardness  and 
of  single  effective  uniform  hardness 

Our  general  analytic  results  predict 
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a)  significant  dispersion/hardness  are  needed  to  achieve  moderate 
levels  of  prompt-effects  survivors  from  attacks  against  population  by  a  few 
thousands  of  megaton-size  weapons 

b)  modest  amounts  of  radiation  shielding  suffice  to  guarantee  that 
no  attack  against  population  can  produce  more  casualties  by  blast  and  fallout 
combined  than  the  casualties  produced  by  blast  alone  from  the  same  number  of 
optimized  airbursts— regardless  of  the  spatial  distribution  of  the  population, 
regardless  of  the  size  and  location  of  the  attack,  regardless  of  the  wind  and 
weather,  regardless  of  the  number  of  casualties 

c)  both  accurate  numerical  values  and  generalized  scaling  laws  to 
understand  trends,  sensiti vites,  regularities. 

Having  indicated  the  origin  of  the  analytic  results  presented  above, 
we  now  turn  to  extending  and  applying  these  results. 

2.2  Effective  Prompt-Effects  Hardness 

We  now  consider  the  realistic  situation  wherein  there  are  M  (mutually 
exclusive  and  collectively  exhaustive)  exposure  environments  for  prompt  effects, 
with  each  such  environment  characterized  by  its  damage  overpressure,  Ap  ,-n  . , 
and  containing  the  fraction  a.  ofthe  total  population,  l<j<M;  we  take  the 

J 

exposure  environments  to  be  ordered  in  terms  of  increasing  damage  overpressure 
so  that  exposure  environments  1  and  M  are  the  "softest"  and  "hardest”,  respec¬ 
tively.  The  damage  distances^  .  is  related  to  the  damage  overpressure 
Ap  cn  h  via  Equations  9  and  10  for  optimized  airbursts  and  surface  bursts, 

S ,DU  > J 


respectively.  Since  our  intent  is  to  explore  the  notion  of  effective  blast 
hardness,  we  restrict  attention  in  this  subsection  to  prompt-effects  casualties. 

By  rank-ordering  cells  of  characteristic  (center-to-corner)  size 
s  and  selecting  aim-points  at  or  above  the  centers  of  the  N  most  densely 
populated  cells,  an  attack  with  N  weapons  produces  prompt-effects  fatalities 
given  by 


where  the  population-weighted  coverage  coefficient  <Q>  is  given  the  obvious 
definition  in  terns  of  the  coverage  coefficient  Q  given  by 


Q(w )  - 


U/2)  w2 

w2  [u  +  arc  cos  u] 


1 


u 


0sw<2 


-1/2 


1  <W<“ 


(31) 


(32) 


It  is  clear  that  Q(w)  is  an  increasing  function  of  w  and  consequently  that 
<Q>  in  Equation  30  is  a  decreasing  function  of  s;  F  in  (3)  is,  of  course,  an 
increasing  function  of  s.  Considered  as  a  function  of  s,  V(s)  has  (as  shown 
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-  1/2 

in  Ref.  1)  for  s  between  Sgg  M  and  2  s^g  ^  a  unique  global  maximum  at,  say, 
s*.  The  effective  damage  distance  and  damage  overpressure  are  given  by 


s50  = 


(33) 


Aps,50  =  Aps^s50^ 


(34) 


In  terms  of  these  effective  values,  we  take  the  prompt-effects  fatalities  to 
be  given  (as  in  Equation  5  above)  by 


n 

n 


fat 

tot 


*  F 


(35) 


In  general,  the  effective  hardness  defined  by  Equations  33  and  34  depends  on 
N  and  F.  We  show  below  that  this  dependence  is  weak.  Moreover,  we  show  that 

A  A 

Equation  35  underestimates  the  "exact"  maximum  V(s=s*)  by  a  negligibly  small 
amount  and  therefore  that  the  combination  of  Equations  33,  34  and  35  provides 
a  valid  representation  of  effective  hardness  and  casualties  for  a  non-uniform 
population. 


Before  investigating  more  complicated  situations,  we  discuss  the 
simplest  case,  where  M  =  1,  i.e.,  there  is  but  a  single  uniform  hardness. 
It  is  clear  from  Equations  33  and  34  that  our  general  approach  specializes 


correctly  in  this  case  to  the  given  single  hardness.  Nevertheless  it  is 
instructive  to  explore  the  variation  of  s*  and  V(s*).  For  concreteness  (and 
without  any  real  loss  of  generality),  in  this  and  all  other  cases  discussed 
in  this  subsection,  we  restrict  attention  to  surface  bursts  and  to  the  distri¬ 
bution  F=Fq  which  is  given  by  Equation  29  and  which  accurately  represents 

American  residential  population.  Our  results  are  summarized  in  Tables  1  and 

2 

2  and  plotted  in  Figures  2  and  3.  We  see  that  as  iNs^/A*  increases,  s*/SgQ 

decreases  steadily  from  1.3502  toward  1  and  that  the  ratio  of  V(s*)  to 

F(tNs,-q/A*)  decreases  steadily  from  1.2289  toward  1.  However,  the  difference 

between  these  two  estimates  of  fatalities  rises  to  a  maximum  of  0.0487  for 
2 

tNs^q/A*  =  0.433  and  then  decreases  steadily.  The  "exact"  maximum  at  s*  greater 
than  s5Q  corresponds,  of  course,  to  targeting  for  an  overpressure  somewhat 
less  than  Ap$  ^  and  accepting  a  smaller  fraction  of  fatalities  per  cell  in 
return  for  targeting  cells  with  a  greater  total  area.  The  approximation 
(underestimation  of  the  "exact"  maximum)  embodied  (via  Equations  33  and  34) 
in  Equation  35  never  significantly  underestimates  the  number  of  fatalities  and 
may  legitimately  be  interpreted  as  enlarging  the  region  in  parameter-space 
wherein  the  approximations  underlying  Equation  35  are  valid. 

To  illustrate  the  more  general  case  of  multiple  exposure  environ¬ 
ments,  we  consider  three  examples  wherein  the  damage  overpressures  Ap  rn  . 
are  the  same  but  the  occupancy  fractions  a.  are  different.  Specifically,  we 

J 

consider  M=4  exposure  environments  with  j  *  1,2, 3, 4  corresponding  to  damage 
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Table  1 


n 

S*/S50 

S50/J* 

<Q>(s=s*) 

.  2 

F(tNs*  /A*) 

V  (s=s*) 

*•# 

10-6 

1.3502 

0 . 7406 

0.8418 

2.422xl0'4 

2.039xl0‘4 

i 

10"5 

1.3501 

0.7407 

0.8419 

1.033xl0"3 

8.696x1  O'4 

s. 

■ » 

1 1 

10"4 

1.3497 

0.7409 

0.8422 

4.397xl0'3 

3.703xl0~3 

.  f 

V*  f, 

10‘3 

1.3481 

0.7418 

0.8435 

0.0186 

0.0157 

l ' 

10‘2 

1.3409 

0.7458 

0.8491 

0.0764 

0.0649 

0.1 

1.3093 

0.7638 

0.8724 

0.2805 

0.2447 

i 

& 

0.2 

1.2865 

0.7773 

0.8881 

0.3925 

0.3485 

X 

0.3 

1.2682 

0.7885 

0.9001 

0.4684 

0.4216 

1 

0.4 

1.2526 

0.7983 

0.9099 

0.5256 

0.4782 

n  ■ 

'¥ 

* 

0.5 

1.2390 

0.8071 

0.9182 

0.5711 

0.5243 

f 

0.6 

1.2268 

0.8152 

0.9253 

0.6085 

0.5630 

* 

0.7 

1.2158 

0.8225 

0.9315 

0.6400 

0.5962 

0.8 

1.2057 

0.8294 

0.9369 

0.6671 

0.6250 

r 

0.9 

1.1966 

0.8357 

0.9418 

0.6906 

0.6504 

‘ 

1.0 

1.1881 

0.8417 

0.9461 

0.7114 

0.6730 

►  , 

*; 

1.5 

1.1539 

0.8666 

0.9623 

0.7870 

0.7573 

2.0 

1.1288 

0.8859  - 

0.9727 

0.8352 

0.8124 

3.0 

1.0942 

0.9139 

0.9847 

0.8933 

0.8796 

r 

4.0 

0.0713 

0.9335 

0.9909 

0.9266 

0.9182 

5.0 

0.0552 

0.9477 

0.9944 

0.9476 

0.9423 

6.0 

1.0434 

0.9584 

0.9965 

0.9617 

0.9583 

;■ 

7.0 

1.0345 

0.9666 

0.9977 

0.9715 

0.9693 

8.0 

1.0277 

0.9730 

0.9985 

0.9784 

0.9770 

9.0 

1.0224 

0.9781 

0.9990 

0.9835 

0.9826 

10.0 

1.0182 

0.9821 

0.9994 

0.9873 

0.9866 
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Table  2 


Comparison  of  Exact  and  Approximate  Values  of  Maximum  Prompt 
Fatalities  -  Single  Exposure  Environment 


tNs50/A* 

F(tNs52/A*) 

V(s=s*) 

_V  (s=s*) 

V  (s=s*) 

F(tNs52/A*) 

10'6 

1.659xl0'4 

2.039xl0~4 

1 .2289 

3.80xl0-5 

10"5 

7.077xl0~4 

8.696xl0"4 

1.2288 

2.61 9x7  0-4 

l(f4 

3.01 5x1 0-3 

3.703xl0“3 

1.2281 

6. 878x1  O'4 

10'3 

0.0128 

0.0157 

1 .2252 

0.0029 

10‘2 

0.0535 

0.0649 

1.2139 

0.0114 

0.1 

0.2090 

0.2447 

1.1710 

0.0357 

0.2 

0.3043 

0.3485 

1.1455 

0.0443 

0.3 

0.3740 

0.4216 

1.1273 

0.0476 

0.4 

0.4296 

0.4782 

1.1132 

0.0486 

0.5 

0.4760 

0.5243 

1.1016 

0.0484 

0.6 

0.5156 

0.5630 

1 .0920 

0.0474 

0.7 

0.5501 

0.5962 

1.0837 

0.0460 

0.8 

0.5806 

0.6250 

1.0765 

0.0444 

0.9 

0.6077 

0.6504 

1.0703 

0.0427 

1.0 

0.6321 

0.6730 

1.0647 

0.0409 

1.5 

0.7250 

0.7573 

1.0445 

0.0323 

2.0 

0.7872 

0.8124 

1.0319 

0.0251 

3.0 

0.8644 

0.8796 

1.0176 

0.0152 

4.0 

0.9088 

0.9182 

1 .0103 

0.0094 

5.0 

0.9365 

0.9423 

1.0063 

0.0059 

6.0 

0.9546 

0.9583 

1.0039 

0.0037- 

7.0 

0.9669 

0.9693 

1.0025 

0.0024 

8.0 

0.9754 

0.9770 

1.0016 

0.0016 

9.0 

0.9815 

0.9826 

1.0010 

0.0010 

10.0 

0.9860 

0.9866 

1.0007 

0.0007 

T  Nsi 


overpressures  of  5,  8,  10  and  35  psi,  respectively,  and  hence  to  DCPA  shelter 
types  G/H/I,  E/F,  B/C/D,  and  A,  respectively.  (For  a  discussion  of  DCPA 
shelter  types,  cf.  Ref.  3.)  For  the  Y  =  1  Mt  surface  bursts  considered  here, 
the  damage  distances  s5Q  j  are  2.869,  2.216,  1.969  and  1.076  mi  for  j  *  1,2, 3,4 
respectively;  here,  and  throughout  this  report,  "mile"  denotes  a  statute  mile 
of  5280  ft  =  1.6093  km.  The  occupancy  fractions  are  selected  from  the  twenty 
combinations  of  four  regions  (Nation,  Northeast,  Midwest,  California)  and  five 
cell  sizes  (length  of  cell  side  corresponding  to  1,  2,  4,  6,  15  min.)  studied 
earlier  (cf.  Table  8,  Ref.  1)..  For  each  of  these  20  combinations,  the  resi¬ 
dential  population  within  each  cell  is  located  in  the  best  (i.e.,  highest 
damage  overpressure)  shelter  available  in  the  cell  so  that,  in  each  cell, 
first  the  35-psi  shelters  are  filled  to  capacity,  then  the  10-psi  shelters, 
etc.;  it  is  assumed  that  everyone  is  located  in  a  shelter  of  at  least  5  psi 
damage  overpressure.  The  three  particular  sets  of  occupancy  fractions  dis¬ 
cussed  here  correspond  to  (California,  1  min.  cell),  (Nation,  2  min.  cell)  } 

» 

and  (Northeast,  15  min.  cell)  and  are  chosen  to  represent  "soft",  "nominal", 
and  "hard"  postures,  respectively.  Although  these  three  examples  are  not 
indicative  of  "softness"  or  "hardness"  in  an  absolute  sense,  they  are  the 
extremes  of  the  20  cases  studied  and  doubtless  do  bracket  the  hardness  of 
national  residential  population  in  "best  available  shelter".  For  descriptive 
simplicity  throughout  this  subsection,  we  continue  to  refer  to  them  as  "soft", 
"nominal"  and  "hard".  The  numerical  characterization  of  the  three  examples 
is  given  in  Table  3.  It  must  be  emphasized  that  the  tabulated  occupancy 
fractions  serve  here  to  characterize  the  three  postures  considered,  and  are 
not  restricted  to  specific  regions  or  cell  sizes. 


Table  3 


Characterization  of  Three  Hardness  Postures 
Multiple  Exposure  Environments 


a.  for 

J 

j 

Aps,50,j  psi 

s50,j  mi. 

Soft 

Posture 

Nominal 

Posture 

Hard 

Posture 

1 

5 

2.869 

0.6884 

0.3471 

0.0073 

2 

8 

2.216 

0.0453 

0.0320 

0.0073 

3 

10 

1 .969 

0.2537 

0.5984 

0.9221 

4 

35 

1.076 

0.0126 

0.0225 

0.0633 

For  these  three  examples,  we  now  present  numerical  results  obtained 

4  2 

by  specializing  Equation  30  for  t=2  and  Equations  1  and  29  for  A*  =  4x10  mi 
corresponding,  respectively,  to  offense-conservative  results  for  the  entire 
nation.  The  behavior  of  V(s)  as  a  function  of  cell  size  is  illustrated  in  Table  4 


Table  4 

V(s)  for  Nominal  Hardness  Posture  and  Several  Values  of  s  and  N 


and  Figure  4  which  list  and  plot,  respectively,  the  magnitude  of  V(s)  as  a 
function  of  s  for  the  nominal  hardness  posture  and  for  several  values  of  N. 

We  see  that,  for  each  values  of  N,  V(s)  rises  to  a  maximum  and  then  decreases 
as  s  increases.  For  each  of  the  three  hardness  postures  and  for  several  values 
of  N,  we  have  determined  the  location,  s*,  and  magnitude,  V(s*),  of  this  maxi¬ 
mum;  the  results  are  displayed  in  Table  5.  Having  determined  s*,  we  invoke 
Equation  33  to  find  s^g  and  F(xNSgg/A*),  with  the  results  shown  in  Table  6. 

We  note  that  as  the  attack  size  increases,  then  s*  decreases  steadily,  s^g/s* 
increases  (i.e.,  s*/s50  decreases)  steadily  toward  1,  and  sgg  decreases  steadily, 
but  relatively  slowly;  the  relative  variation  in  s^g  is  three  to  five  times 

A 

less  than  that  of  s*.  As  N  increases,  the  difference  between  the  "exact"  maxi- 
mum  V (s* )  and  F(xNs50  /A*)  increases  to  a  maximum  and  then  decreases  steadily; 
the  maximum  value  of  this  difference  is  relatively  small  (less  than  five  per¬ 
cent).  For  the  most  part,  this  behavior  is  precisely  that  demonstrated  earlier 
for  the  case  of  single  uniform  hardness.  These  results  are  depicted  in  Figures 
5,  7,  and  9  which  present  s*  and  s5g  versus  N  for  the  soft,  nominal,  and  hard 
postures,  respectively,  and  in  Figures  6,  8  and  10  which  present  V(s*)  and 
F(tNS5q/A*)  in  like  fashion. 

Having  determined  the  damage  distance  s^g,  we  obtain  the  corresponding 
damage  overpressure  from  Equation  34.  In  Table  7,  we  summarize  the  results  for 
Sjjq  and  aps  gg  for  all  three  hardness  postures  and  for  several  values  of  N; 
these  results  for  s^g  and  Ap$  ^g  are  plotted  in  Figures  11  and  12,  respectively. 
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Table  5 


Location  and  Magnitude  of  V(s* )  for  Three  Hardness  Postures 


Posture 

!  N 

t  - 

<Q>(s=s*) 

V (s=s*) 

Soft 

30 

3.731 

0.7426 

0.0808 

102 

3.679 

0.7549 

0.1619 

0.1222 

3xl02 

3.582 

0.7777 

0.2891 

0.2248 

103 

3.372 

0.8244 

0.4909 

msm 

3.085 

0.8828 

0.7006 

0.6185 

104 

2.704 

0.9425 

0.8870 

0.8360 

Nominal 

30 

3.119 

0.7469 

0.0650 

0.0485 

102 

3.022 

0.7757 

0.1288 

0.0999 

2.784 

0.8517 

0.1873 

103 

2.759 

0.8598 

0.4080 

0.3508 

■a 

2.648 

0.8895 

0.6302 

0.5606 

104 

2.375 

0.9466 

0.8431 

0.7981 

Hard 

30 

2.647 

0.8114 

0.0532 

102 

2.632 

0.81 70 

0.1094 

0.0893 

KUH 

2.603 

0.8275 

0.2040 

0.1688 

io3 

2.536 

0.8504 

0.3760 

0.3197 

3xl03 

2.419 

0.8866 

0.5883 

0.5216 

a 


0.8178 


Table  6 

Effective  Damage  Distance  and  Casualties  for  Three  Hardness  Postures  \ 

Cj 


Posture 

N 

s*,mi . 

s50/S* 

s50’mi  • 

i  F(tNs|0/A*) 

I  ...  _! 

Soft 

30 

3.731 

0.6876 

2.566 

0.0512 

C 

102 

3.679 

0.6932 

2.551 

0.1054 

3xl02 

3.582 

0.7036 

2.520 

0.1967 

103 

3.372 

0.7290 

2.458 

0.3645 

3xl03 

3.085 

0.7726 

2.384 

0.5816 

r 

104 

2.704 

0.8367 

2.262 

0.8248 

Nominal 

30 

3.119 

0.6896 

2.151 

0.0412 

102 

3.022 

0.7027 

2.123 

0.0846 

C 

3xl02 

2.784 

0.7477 

2.082 

0.1582 

103 

2.759 

0.7538 

2.079 

0.3073 

3xl03 

2.648 

0.7786 

2.062 

0.5161 

104 

2.375 

0.8423 

2.001 

0.7750 

(_ 

Hard 

30 

2.647 

0.7211 

1.909 

0.0355 

102 

2.632 

0.7244 

1.906 

0.0742 

3x1 02 

2.603 

0.7310 

1.903 

0.1425 

(_ 

103 

2.536 

0.7467 

1.894 

0.2785 

3xl03 

2.419 

0.7760 

1.877 

■  0.4752 

104 

2.222 

0.8270 

1.837 

0.7382 

< 


gure  5.  s*  and  s 


Table  7 


Effective  Damage  Distance  and  Overpressure  for  Three  Hardness  Postures 


Sjjg.mi.  for 

Aps,50’  pSi 

for 

Soft 

Posture 

Nomi na 1 
Posture 

Hard 

Posture 

Soft 

Posture 

Nomi nal 
Posture 

Hard 

Posture 

2.566 

2.151 

1.909 

6.11 

8.46 

10.63 

2.551 

2.123 

1.906 

6.18 

8.67 

10.65 

2.520 

2.082 

1.903 

6.31 

9.00 

10.69 

2.458 

2.079 

1.894 

6.60 

9.02 

10.79 

2.384 

2.062 

1.877 

6.99 

9.16 

10.97 

2.262 

2.001 

1.837 

7.70 

9.71 

11.44 

Since  Aps  ^  s”n  with  n  *  2  for  this  range  of  overpressures,  the  variation  of 
damage  overpressure  is  approximately  twice  that  of  the  damage  distance.  How¬ 
ever,  the  variation  of  both  is  relatively  small.  For  distributions  which 
are  flatter  than  the  distribution  F  =  F0  employed  here,  the  variation  is  even 
less.  We  have  shown  (cf.  Ref.  1)  that  a  defense-conservative  upper  bound  on 
SgQ  is  obtained  by  taking 


S50  =  s50 


= $5°'0 


Substitution  of  values  from  Table  3  into  Equation  36  gives  s^q  s  2.625  mi., 
2.316  mi.  and  1.936  mi.  for  the  soft,  nominal,  and  hard  postures  considered 
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gure  11. 


here,  with  corresponding  values  for  Aps  ^  of  5.9  psi ,  7.4  psi,  and  10.3  psi , 
respectively.  Comparison  of  these  values  with  those  tabulated  in  Table  7 
shows  that  Equation  36  provides  a  rather  good  approximation  to  s^g  and  hence 
to  Aps  For  many  applications  requiring  a  defense-conservative  estimate, 
the  simplicity  of  Equation  36  in  determining  s^g  and  Ap$  ^g  outweighs  the 
loss  in  accuracy. 

In  sum,  we  have  shown  above  that  our  prescription  (embodied  in 
Equations  33,  34,  35)  for  prompt-effects  hardness  and  casualties  produces 
stable,  relatively  invariant  values  of  effective  hardness  and  approximates 
closely  the  "exact"  maximum  number  of  casualties  achievable. 

2.3  Parametric  Characterization  of  Relocation 

It  is  clear,  from  both  common  sense  and  from  the  precise  results  of 
Equation  23,  that  passive  defense  measures  are  perforce  combinations  of  the 
generic  techniques  of  hardening  and  dispersion.  We  have  discussed  above  the 
characterization  of  hardness,  and  now  turn  to  considering  dispersion.  Special 
interest  attaches  to  relocation  schemes  wherein  people  are  evacuated  to  loca¬ 
tions  where  at  least  some  people  already  reside.  In  this  subsection  we  pre¬ 
sent  a  parametric  characterization  of  such  relocation  policies.  We  first 
present  results  for  relocation  within  a  single  region  and  then  in  view  of 
these  results,  consider  nationwide  relocation  with  "self-contained"  relocation 
of  individual  regions. 
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Our  discussion  of  relocation  is  couched  in  terms  of  the  distributions 
(rank-ordered  by  population  density  and  presenting  cumulative  population  versus 


cumulative  area)  defined  by  the  spatial  distributions  of  population  before 
and  after  relocation.  Unless  otherwise  noted,  we  take  the  r3nk-ordered  distri¬ 
bution  of  residential  (i.e.,  pre-relocation)  population  in  each  region  to  be 
given  by  Equations  1  and  29  (i.e.,  by  F  =  Fq  with  x  =  0.63)  with  Aqcc  =  10A*; 
thus  each  region  is  characterized  by  its  total  population,  ^tot>  and  charac¬ 
teristic  area,  A*.  Since  Fo()0)  =  0.986,  use  of  Fq  with  Aqcc  =  10  A*  implies 
that  1.4  percent  of  the  population  lives  "outside  the  populated  area".  Al¬ 
though  this  fact  is  clearly  not  truly  meaningful,  it  suggests  that  some  caution 
must  be  exercised  in  defining  relocated  populations  in  terms  of  F  .  We  adopt 
the  convention  that  only  the  10A*  region  containing  Fq( 10)  =  98.6  percent  of 
the  total  population  is  a  candidate  for  receiving  evacuees.  Moreover,  we 
recognize,  within  the  populated  area,  three  regions  corresponding  to  decreasing 
ranges  of  population  density  for  residential  population:  evacuation  region, 
from  which  at  least  some,  and  generally  most,  people  are  relocated;  buffer 
region,  neither  into  nor  out  of  which  relocation  occurs;  potential -host 
region,  i.e.,  the  populated  area  remainder  which  is  a  candidate  for  receiving 
evacuees.  We  distinguish  between  the  potential -host  region,  so  defined,  and 
the  host  region,  which  is  the  region  which  actually  receives  evacuees;  if 
the  two  regions  coincide,  then  the  hosting  is  said  to  be  complete. 

Relocation  reduces  vulnerability  insofar  as  the  rank-ordered  distri¬ 
bution  corresponding  to  the  relocated  population  is  flatter  than  the 
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original  distribution.  The  population  density  is  lessened  in  the  evacuation 
region,  is  unchanged  in  the  buffer  region,  and  necessarily  is  increased  in 
the  host  region.  Hence  the  key  to  achieving  a  relatively  flat  distribution 
for  the  relocated  population  is  control  of  the  post-relocation  population 
density  within  the  host  region.  We  therefore  consider  relocation  schemes 
wherein  the  population  density  in  the  host  region  is  constrained  not  to  exceed 
a  prescribed  value,  w*.  A  further,  pragmatic  restriction  arises  from  the 
finite  capacity  of  the  host  region  to  accomodate  evacuees;  an  index  of  this 
restriction  is  the  hosting-ratio  K  defined,  for  any  portion  of  the  host 
region,  by 


r,  _  population  after  receiving  evacuees 
population  before  receiving  evacuees 


We  consider  relocation  schemes  wherein  the  hosting-ratio  in  the  host  region 
is  constrained  not  to  exceed  a  prescribed  value,  K*. 

We  now  summarize  the  resulting  definition  of  relocation  schemes 
which  evacuate  people  to  locations  where  at  least  some  people  already  reside: 

i)  the  evacuation  region  consists  of  the  most  densely  populated  area  of  A, 
and  the  fraction  y  of  the  population  of  the  evacuation  region  is  relocated; 

ii)  the  buffer  region  consists  of  the  locations  where  the  residential  popu¬ 
lation  density  is  less  than  in  the  evacuation  region,  but  greater  than  u*; 
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iii)  within  the  potential-host  region,  the  locations  with  the  highest 
residential  population  density  receive  as  many  evacuees  as  possible,  subject 
to  the  constraints  that  the  population  density  after  relocation  not  exceed 
to*  and  that  the  hosting  ratio  not  exceed  K*.  The  buffer  region  is  repre¬ 
sented  before  relocation  by  the  portion  of  F  *  Fq  extending  from  A  =  A  to 
A  =  A  where,  in  view  of  Equation  2, 


ntot 

A*  h  yA* j 


=  far 


(38) 


Likewise,  the  potential -host  region  is  represented  before  relocation  by  two 

A  _ 

^  'V  % 

portions  of  F  =  Fq  extending  from  A  =  A  to  A  and  from  A  =  A  to  A  =  10A*, 
respectively,  where  A  is  given  by 


n 


tot 


A* 


(39) 


Clearly,  the  portion  extending 'from  A  to  A  is  limited,  after  relocation,  by 

the  constraint  that  the  population  density  not  exceed  u>*,  while  the  portion 

% 

extending  from  A  to  10A*  is  limited  by  the  constraint  that  the  hosting-ratio 

not  exceed  K*.  (Ref.  1  illustrates  this  generic  relocation  scheme  in  a  single 

specific  example  characterized  by  takingthe  region  to  be  the  entire  nation, 

'  3  2 

A  =  1.0A*,  y  =  0.8,  u*  *  10  people/mi  ,  K*  =  7;  the  resulting  relocated 

distribution  is  denoted  F  .) 

con 
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For  a  given  co*  and  K*,  the  potential -host  region  has  a  maximum 


hosting  capacity  (i.e.,  excess  of  population  which  can  be  contained  over  and 
above  the  residential  population).  This  maximum  capacity  occurs  for  complete 
hosting  and  is  given  by 


[f(10)-f  (fcj  ] 


capacity  =  oi*[A-A]  +  K*ntQt  |F(10)-F 


‘  ntot[F(1°)-F(^r 


(40) 


In  dimensionless  form  we  have  this  capacity  (normalized  by  total  population) 
as  r  where,  from  Equations  38,  39,  40 


r ' F' (^)[^'Mtii*[F(,0)'F  (**)] 

-[f(10)-f(^)]  (41) 

In  view  of  Equations  38  and  39,  we  also  have 

r  •  {  f'(tf)  [f  -  f]+  (jpr)}  -  [F(10>-F  [wj]  <42> 

The  hosting  demand  (i.e.,  the  additional  population  to  be  relocated  within 
the  host  region)  is,  of  course,  simply  the  number  of  evacuees.  In  dimensionless 


terms,  the  demand  (normalized  by  total  population)  is  denoted  S  and  given 
by 


S  =  y  F(jj*)  (43) 

In  order  for  the  relocation  to  be  feasible,  the  hosting  capacity  must  exceed 
the  hosting  demand,  i.e., 


r  >  S  (44) 

We  now  investigate  the  conditions  for  which  this  fundamental  in¬ 
equality  is  satisfied.  We  restrict  attention  to  relocations  for  which  there 
is  no  buffer  region,  so  that  A  =  A;  this  serves  both  to  reduce  the  already 
large  number  of  parameters  and  to  create  the  maximum  possible  hosting  demand. 
(The  results  presented  below  are  readily  extended  to  include  a  finite  buffer 
region;  specific  buffer  regions  considered  in  the  past  are  difficult  to 
characterize  in  intrinsic  terms,  but  generally  have  relatively  small  total 
population.)  Our  approach  is  to  characterize  relocations  with  complete 
hosting;  the  resulting  post-relocation  rank-ordered  distributions  are  denoted 
by  G  (to  avoid  confusion  with  the  pre-relocation  distribution  F  =  F  ). 

We  proceed  by  first  determining  the  minimum  permissible  u*  for 


given  K*  and  y.  Determining  this  minimum  permissible  u;*  is  equivalent  (via 


1 


Equation  38)  to  determining  the  maximum  permissible  A  and  hence  is  equivalent 
(via  Equation  39)  to  determining  the  maximum  permissible  A.  Of  the  various 
alternative  ways  to  achieve  these  equivalent  ends,  perhaps  the  simplest  is 
interpreted  as  an  iterative  procedure  (with  convergence  guaranteed  by  the 
fact  that  the  procedure  is  demonstrably  equivalent  to  inverting  a  monotonic 
function)  wherein  one  i)  picks  a  value  of  A,  ii)  thereby  determines  A  since 
Equations  38  and  39  imply  that 


F'(;pr)  ■  K*  F(jjr) 


(45) 


iii)  thereby  determines  r  from  Equation  41  or  42,  iv)  compares  r  with  S  deter- 

A  A 

mined  from  Equation  43  with  A  =  A,  v)  iterates  to  obtain  r  =  S.  We  denote 

A  ^ 

the  results  in  terms  of  the  parameters  a  and  a  which  give  the  maximum  per¬ 
il  'v# 

missible  values  of  A  and  A  via 


2  _  maximum  permissible  A 

a  “jpr - 


(46) 


a  * 


% 

maximum  permissible  A 
A* 


(47) 


It  obviously  follows  that 


minimum  -permissible  u*  = 


‘tot 


A* 


F’(a)  =  K* 


tot 

A* 


F'fl) 


(48) 
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The  results  are  summarized  in  Tables  8  and  9  and  Figures  13,  14, 
and  15.  Tables  8  and  9  list  critical  values  of  hosting  parameters  (e.g., 

A 

A  ^  A 

a,  a,  and  functions  thereof)  as  functions  of  K*  for  y  =  1.0  and  0.8,  i.e., 
for  100  percent  and  80  percent  evacuation,  respectively.  Figures  13  plots 
the  minimum  permissible  hosting  population  density  m*  as  a  function  of  K* 
for  y  =  1.0  and  0.8;  alternatively.  Figure  13  depicts  the  boundary  in  popu¬ 
lation-density/hosting-ratio  space  between  feasible  and  infeasible  relocations. 

In  like  fashion.  Figures  14  and  15  present  the  variation  with  K*  of  the  maxi¬ 
mum  evacuated  population,  yF(a),  and  maximum  evacuation  area,  a,  respectively; 
these  Figures  may  also  be  interpreted  as  defining  boundaries  in  the  appropriate 
parameter  space.  We  note  the  results  that,  for  K*  =  7,  then  the  minimum  per¬ 
missible  hosting  population  density  w*  is  0.1741  ntot/A*  and  0.1433 
for  y=  1.0  and  y=  0.8,  respectively;  the  analogous  values  of  maximum  evacuated 
population  are  0.6945  ntQt  and  0.5860  while  the  analogous  values  of 

maximum  evacuated  area  are  1.3106  A*  and  1.5511  A*,  respectively.  We  further 

note  the  consequence  that  for  y  -  0.8  and  K*  =  7,  a  hosting  population  density 
3  2 

w*  =  10  people/mi  implies  a  feasible  relocation  if,  and  only  if,  &tot/A*  is 
less  than  (0.1433)"1  10^  people/mi2  =  6.98x10^  people/mi2  =  7x10^  people/mi2. 

It  is  apparent  from  Figures  13  through  15  that  for  hosting-ratios  larger 
than  6  to  8,  the  critical  values  of  population  density,  evacuated  population  and 
evacuation  area  tend  to  flatten  out;  thus  the  dependence  upon  hosting-ratio 
is  weak,  once  the  hosting-ratio  reaches  approximately  7.  Additional  insight 
into  the  behavior  of  the  relocation  is  obtained  by  examining  how  the  population 
after  relocation  is  distributed  among  the  three  (mutually  exclusive  and 


Table  8 


Table  9 

Critical  Values  of  Hosting  Parameters  for  80%  Evacuation 


Figure  14.  Maximum  Permissible  Evacuated  Population  vs  Hosting-Ratio 


Maximum  Permissible  Evacuation  Area  vs  Hostinq-Ratio 


collectively  exhaustive)  regions:  evacuation  region,  host-region  with  constant 
population  density,  host-region  with  constant  hosting-ratio.  Figures  16  and 
17  present  the  population  of  these  three  regions  as  a  function  of  hosting-ratio 
for  y  =  1.0  and  y  =  0.8,  respectively,  for  y  =  1.0,  the  post-relocation  popula¬ 
tion  of  the  evacuation  region  is,  of  course,  zero.  We  see  from  these  Figures 
that  the  post-relocation  population  of  the  evacuation  region  depends  weakly 
(e.g.,  varying  between  10  and  16  percent  for  y  =  0.8)  on  hosting-ratio,  while 
the  post-relocation  host  region  population  with  constant  population  density  in¬ 
creases  steadily  with  increasing  hosting-ratio;  necessarily,  the  remaining  pop¬ 
ulation  (i.e.,  host  region  with  constant  hosting-ratio)  decreases  steadily  with 
increasing  hosting  ratio.  Although  these  breakdowns  of  the  total  post-relocation 
population  are  instructive,  it  is  even  more  informative  to  determine  the  rank- 
ordered  (by  population  density)  distributions  G  which  are  defined  by  the  re¬ 
located  population. 

It  is  clear  that,  for  complete  relocation  with  r  =  S,  the  rank-ordered 
distribution  G  begins  with  that  portion  of  the  evacuation  region  which  even  after 
relocation  has  population  density  greater  than  w*.  The  next  portion  of  G  has 
uniform  population  density  u*.  In  rank-ordering  the  balance  of  the  relocated 
population  (i.e.,  the  remainder  of  the  evacuation  region,  plus  the  host  region 
with  constant  hosting-ratio)  it  is  necessary  to  recognize  that  the  smallest 
population  density  in  the  evacuation  region  is  (I-y)cd*,  while  the  smallest  pop¬ 
ulation  density  in  the  host  region  is  K*F' (1 0)ntQt/A*,  where  F ’ (10)  =  3 . 773x1 0~3 . 


Breakdown  of  Post-Relocation  Population  vs  Hostinrj-Ratio  for  100%  Evacuation 


Figure  17.  Breakdown  of  Post-Relocation  Population  vs  Hosting-Ratio  for  80'X  Evacuation 


For  a  given  y,  the  result  is  that  for  sufficiently  small  hosting-ratio  K*, 
the  least  densely  populated  cells  are  in  the  host  region,  i.e.,  K*F‘ (lG)fitot/A* 

<  (1-y)w*;  for  sufficiently  large  hosting-ratio,  the  converse  holds.  As 
may  be  seen  from  Table  9  for  y  -  0.8,  the  transition  occurs  for  K*  between  7 
and  8.  Thus  for  80  percent  evacuation,  the  least  densely  populated  portion  of  the 
relocated  population  is  located  in  the  host  region  for  K*  <  7  and  in  the  evacua¬ 
tion  region  for  K*  >  8.  In  general,  then,  after  the  portion  with  constant  pop¬ 
ulation  density  u*,  the  rank-ordered  distribution  6  has  a  portion  containing 
people  from  the  evacuation  region  and  from  the  host  region  with  population 
density  between  u>*  and  the  greater  of  the  two  minimum  population  densities  of 
these  two  regions.  The  final  portion  of  G  then  contains  the  balance  of  the 
residents  of  the  region  with  smaller  minimum  population  density.  The  precise, 
quantitative  characterization  of  the  distribution  G  is  presented  in  Appendix  A, 
which  may  be  consulted  for  details  and  numerical  examples. 

The  preceding  discussion  has  considered  a  single  region  exclusively 

and  the  relocation  of  the  region's  population  within  the  region  itself.  Such 

regionally  self-contained  relocations  are  an  integral  part  of  relocation  on  a 

larger  (e.g.,  national)  scale.  Therefore,  the  results  presented  above  have 

significant  consequences  for  larger  scale  relocation.  For  example,  we  have  seen 

above  that  for  80  percent  evacuation  and  a  critical  hosting-ratio  of  7  then  a 
.  2 

critical  population  density  of  1,000  people/mi  implies  a  feasible  relocation 
if  and  only  if  the  region's  ^tot/A*  is  less  than  7.0xl03  people/mi2.  Con¬ 
sequently,  (frpm  the  demographic  statistics  of  Table  3,  Ref.  1)  we  see 
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that  relocation  with  these  values  of  the  hosting  parameters  is  feasible  for  the 

entire  CONUS  and  for  the  Midwest  (with  nt()t/A*  of  5.7xl03  and  6.8xl03  people/mi2, 

respectively),  but  not  for  the  Northeast  and  California  (with  q  t/A*  of 
3  3  2 

9.8x10  and  10.2x10  people/mi  ,  respectively).  Thus  the  specific  relocation 

described  in  Ref.  1  is  feasible  on  a  national  scale,  but  has  the  drawback  that 

residents  of  the  Northeast  and  California  must  be  relocated  outside  their  own 

regions.  The  definition  of  this  specific  relocation  has  been  repeated  above 

and  can  be  paraphrased  as  calling  for  i)  80  percent  evacuation  of  the  locations 

3  2 

wherein  the  residential  population  density  exceeds  1.2x10  people/mi  ,  ii) 

a  buffer  region  consisting  of  residential  population  with  population  density 

3  2  3  2 

between  1.2x10  people/mi  and  u*  =  1.0x10  people/mi  ,  iii)  a  potential  host 

region  wherein  the  population  density  after  relocation  does  not  exceed  u.* 

A 

and  the  hosting-ratio  does  not  exceed  K*  =  7. 

For  the  Northeast  (with  =  9.8xl03  people/mi2),  the  values  of 

3  2  « 

u*  =  10  people/mi  and  K*  =  7  imply  that  a  =  2.028  and  a  =  6.039;  substitution 

of  these  values  into  Equations  41  or  42  gives  a  maximum  capacity  of  r  =  0.4293. 

However,  the  evacuation  of  residential  locations  with  population  density  greater 
3  2 

than  1.2x10  people/mi  implies  that  A  is  1.792  times  the  characteristic  area  A* 
for  the  Northeast;  hence,  from  Equation  43,  the  demand  (i.e.,  the  number  of 
evacuees)  is  0.6112.  Consequently,  the  demand  exceeds  the  capacity  by  (0.6112- 
0.4293)  =  0.1819  times  the  regional  population  of  54.9  million  people;  thus 
10.0  million  residents  of  the  Northeast  cannot  be  relocated  within  the  Northeast. 
Moreover,  the  fraction  (0.1819/0.6112)  =  29.8  percent  of  the  evacuees  can't  be 
relocated  within  the  Northeast. 


For  California,  the  analogous  calculations  give  the  results  that  a  =  2.090, 
a  =  6.143,  r  =  0.4127  while  S  =  0.6167  for  A  equal  to  1.850  times  the  regional 
characteristic  area  A*.  Consequently,  demand  exceeds  capacity  by  (0.6167- 
0.4127)  =  0.2040  times  the  regional  population  of  21.1  million;  thus  4.3 
million  California  residents  cannot  be  relocated  within  California.  Likewise, 
the  fraction  (0.2040/0.6167)  =  33.1  percent  of  the  evacuees  cannot  be  relocated 
within  California.  In  sum,  then,  we  see  that  for  both  the  Northeast  and 
California,  some  one-fifth  of  the  total  regional  population  and  some  one- 
third  of  the  total  number  of  evacuees  cannot  be  relocated  within  their  own 
region. 

We  now  illustrate  how  the  abstract  results  derived  above  can  be 
applied  to  define  a  nationwide  relocation  for  which  relocation  of  the  Northeast 
and  California  is  self-contained.  We  again  consider  relocation  with  evacua¬ 
tion  fraction  y  =  0.8  and  critical  hosting-ratio  K*  of  7;  as  noted  above, 
hosting-ratios  larger  than  7  have  limited  payoff  and  cause  the  entire  host 
region  to  be  more  densely  populated  than  some  portions  of  the  evacuation  region. 
Since  the  Northeast  and  California  have  very  nearly  the  same  residential  popu- 
lation  density  ft^/A*  people/mi  ),  we  treat  them  as  a  single  region 

with  total  population  of  76  million  and  characteristic  area  A*  of  7.5xlC  mi  , 

and  we  consider  the  remainder  of  the  nation  as  a  second  region  with  total  pop- 

3  2 

ulation  of  136  million  and  characteristic  area  A*  of  30.6x10  mi  ;  to  within 
roundoff,  these  values  preserve  the  total  population  and  total  occupied  areas 
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of  the  regions  (given  in  Ref.  1,  Table  3,  based  on  DCPA  extrapolation  of  1970 

census  data),  while  taking  AQCC  =  10A*.  These  values  imply  that  is 

2  2 

10,133  people/mi  and  4,444  people/mi  for  the  two  regions.  By  demanding  that 

the  relocation  of  the  artificial  composite  region  consisting  of  the  Northeast 

plus  California  be  self-contained,  we  obviously  guarantee  that  the  relocation 

of  the  actual,  geographically  distinct  regions  be  individually  self-contained. 

From  Table  9  (or  Figure  13),  we  find  that,  for  y  =  0.8  and  K*  =  7,  the  minimum 

permissible  hosting  population  density  for  self-contained  relocation  is  0.1433 

(10,133  people/mi  )  =  1452  people/mi  for  the  Northeast  and  California.  Con- 

2 

sequently,  a  national  relocation  with  w*  =  1452  people/mi  and  K*  =  7  will 
permit  residents  of  the  Northeast  and  California  to  be  relocated  in  their  own 
regions,  with  complete  hosting. 

The  national  relocation  thus  inferred  calls  for  i)  evacuating  80 

percent  of  the  population  in  all  locations  where  the  residential  population 

2 

density  exceeds  1452  people/mi  ,  ii)  having  no  buffer  region,  iii)  restricting 

2 

the  population  density  and  hosting-ratio  not  to  exceed  1452  people/mi  and  7, 
respectively,  in  the  host  region.  By  construction,  the  hosting  of  the  North¬ 
east  and  California  is  complete,  with  parameters  given  by  Table  9.  For  the 
region  consisting  of  the  balance  of  the  nation,  we  find  that  the  evacuation 
region  has  a  residential  population  of  0.5476  times  the  total  regional  popu¬ 
lation  and  an  area  of  0.6919  times  the  region's  characteristic  area.  Consequently, 
the  number  of  evacuees  is  0.8  (0.5476)  =  0.4380  times  the  regional  population. 
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We  find  that  these  evacuees  are  accomodated  in  the  host  region  with  uniform 

? 

population  density  of  1452  people/mi  over  an  area  of  2.3107  times  the 
characteristic  area  of  the  region.  In  the  composite  rank-ordered  distribu¬ 
tion  for  the  entire  (relocated)  nation,  we  find  that  4.20  percent  of  the 

o 

national  population  has  population  density  greater  than  1452  people/mi  and 

3  2 

covers  an  area  of  3.22x10  mi  ,  while  67.07  percent  of  the  national  popula- 

2 

tion  has  population  density  equal  to  1452  people/mi  and  covers  an  area  of 
3  2 

97.92x10  mi  .  The  balance  of  the  population  (less  than  30  percent  of  the 

total)  is  contained  in  the  evacuation  region,  in  the  constant  hosting-ratio 

portion  of  the  host  region  for  the  Northeast  and  California,  and  in  the 

remainder  of  the  potential  host  region  for  the  rest  of  the  nation;  the  last 

of  these  are  at  residential  population  density.  We  call  this  rank-ordered 

distribution  Fcr  and  plot  it  in  Figure  18  along  with  the  distribution  F 

ic  con 

derived  earlier  (Ref.  1);  in  plotting  Figure  18,  we  take  the  characteristic 

,  .  4  2 

area  A*  for  the  nation  to  be  4x10  mi  ,  i.e.,  the  value  used  earlier  in  con¬ 
structing  Fcon-  We  also  include  in  Figure  18  the  distribution  F  =  Fq  corres¬ 
ponding  to  residential  population  and  the  absolutely  limiting  distribution  F  =  F^ 
corresponding  to  uniform  population  density  over  the  entire  populated  area. 

It  is  apparent  from  Figure  18  that  requiring  the  relocation  of  the 
Northeast  and  California  to  be  self-contained  has  steepened  the  rank-ordered 
distribution  of  the  relocated  population.  Moreover,  by  comparison  with  the 
irreducible  limit  of  F  =  F^  with  its  uniform  population  density  of  approximately 
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500  people/mi  ,  the  distributions  FcQn  and  F^  have  the  bulk  (greater  than 

70  percent  of  the  national  population)  at  population  densities  of  approximately 

2 

1,000  and  1,500  people/mi  ,  respectively,  i.e.,  approximately  2  and  3  times 
the  irreducible  limit.  Despite  the  idealizations  in  our  derivation  of  the 

distributions  F_r  and  F _ ,  it  is  highly  plausible  that  no  rank-ordered  nation- 

wide  distribution  flatter  than  something  intermediate  between  F<-c  and  FcQn  can 
be  produced  by  any  realistic  relocation  scheme  which  evacuates  people  to  loca¬ 
tions  where  at  least  some  people  already  reside  and  which  guarantees  self- 
contained  relocation  for  regions  of  reasonable  size.  Moreover,  the  detailed 
definition  of  any  such  relocation  scheme  must  recognize  the  inherent  capacity 
vs.  demand  limitations  explored  above  and  will  be  simplified  by  applying  the 
specific  abstract  results  derived  above. 


2.4  Applications  and  Examples 

Section  2.1  above  has  summarized  a  number  of  fundamental  results 
concerning  the  structure  of  casualty  prediction  problems.  Because  of  the 
analytic  simplicity  of  these  results,  it  is  straightforward  to  characterize 
parametrically  the  interrelationships  among  the  problem's  ingredients.  In 
this  subsection,  rather  than  constructing  an  encyclopedic  catalog  of  such 
parametric  results,  we  present  numerical  examples  and  analytic  extensions 
choosen  to  illustrate  some  basic  properties  of  passive  defense. 
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We  begin  by  exploring  the  relative  influence  of  hardening  and 
dispersion.  We  consider  prompt-effects  casualties  exclusively;  equivalently, 
we  consider  combined-effects  casualties  for  populations  whose  weighted  pro¬ 
tection  factor  is  large  enough  to  satisfy  the  condition  of  Equation  25.  We 
consider  airbursts  exclusively,  optimized  for  the  prompt-effects  hardness 
aps  go*  We  have  seen  above  that  the  prompt-effects  hardness  for  the  residential 
population,  even  in  "best  available  shelter",  lies  between  5  and  10  psi;  we 
typically  take  7  psi  as  a  representative  value  for  the  prompt-effects  hardness. 
The  spatial  arrangement  of  the  population  is  represented  by  the  distribution  F, 
rank-ordered  by  population  density;  we  have  seen  above  that  F  =  Fq  for  the 
residential  population,  while  the  most  defense-optimistic  F  lies  between  Fsc 
and  Fcon  for  population  relocated  to  places  where  at  least  some  people  already 
reside. 

Figure  19  (based  on  Equation  23)  presents,  as  a  function  of  blast- 
2/3 

equivalent  megatons  (i.e.,  NY  ),  defense-conservative  estimates  of  the  prompt- 
effects  hardness  needed  to  achieve  survivor  fractions  of  1/2  and  2/3  for  the 
residential  and  relocated  (according  to  F-r  and  F  )  populations.  It  is 
apparent  that  even  for  the  relocated  populations,  significant  hardness  is 
required  for  not  unreasonable  attack  sizes.  An  alternative  viewpoint  is  pre¬ 
sented  in  Figure  20  (based  on  Equation  5)  which  plots  a  defense  conservative 
estimate  of  fatalities  as  a  function  of  blast-equivalent  megatons  for  a  pop¬ 
ulation  with  aps  5q  =  7  psi  and  with  spatial  distribution  corresponding  to 


ze  for  7-psi  Population,  Residential  and  Relocated 


of  Blast-Equivalent  Megatons 


Fq  and  f:con*  From  one  viewpoint  (i.e.,  read  "vertically"),  Figure  20  indicates 
that  for  a  given  attack  size,  relocation  can  reduce  fatalities  appreciably; 
from  another  viewpoint  (i.e.,  read  "horizontally").  Figure  20  indicates  that 
for  a  given  fatality  level,  relocation  can  increase  the  required  attack  size 
by  no  more  than  a  few  thousand  megatons.  In  Figure  21  (based  on  Equation  5), 
we  illustrate  the  effect  of  increasing  hardness  without  relocation,  by  plotting 
fatalities  versus  blast-equivalent  megatons  for  prompt-effects  hardnesses  of 
35  and  50  psi;  Figure  21  presents  both  offense-conservative  and  defense- 
conservative  estimates  so  that  the  relatively  small  difference  between  the  two 
bounds  may  be  recognized.  The  inevitable  conclusion  from  these  results  is 
that  to  achieve  moderate  survivor  fractions  for  attacks  against  population  with 
several  thousand  blast-equivalent  megatons,  significant  hardening  is  required, 
even  with  the  most  defense-optimistic  assumptions  about  relocations  which 
evacuate  people  to  locations  where  at  least  some  people  already  reside. 

We  now  turn  to  exploring  the  influence  of  the  delayed-effects  hard¬ 
ness,  <p>.  For  combined  (i.e.,  prompt-plus-delayed  effects)  casualties  to 
exceed  the  (exclusively  prompt-effects)  casualties  from  an  all-airburst  attack, 
it  is  necessary--but  by  no  means  sufficient--that  <p>  be  sufficiently  small 
that  the  condition  of  Equation  25  is  not  satisfied.  Even  for  such  small  values 
of  <p>,  the  excess  A  of  combined  casualties  (from  a  mix  of  surface  and  air- 
bursts)  over  prompt  casualties  (from  a  pure  airburst  attack)  is  bounded  via 
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tot 


Figure  21 . 


Fatalities  vs  Attack  Size  for  Residential  Roculation, 
35  and  50  psi  Hardness 
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The  value  of  delayed-effects  hardness  for  which  Equation  25  holds  as  an  equality 

may  naturally  be  called  the  breakeven  value  of  protection  factor,  since  this 

value  gives  <j>  =  0  in  Equation  49,  while  larger  values  give  <j>  <  0  and  smaller 

values  give  <j>  >  0.  In  our  discussion  below,  we  restrict  attention  to  values 

of  <p>  less  than  the  breakeven  value  and  to  F  =  F  so  that  the  results  are 

o 

applicable  to  residential  population. 

We  begin  by  listing  in  Tables  10,  11,  12  (for  Ap$  5Q  =  5,  7,  10  psi, 
respectively)  the  values  of  the  bound  given  by  Equation  49  for  Y  =  1  Mt  and 
various  attack  sizes  and  values  of  protection  factor.  In  interpreting  these 
tables,  it  is  instructive  to  recall  that 
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for  A*  =  4x10  mi  and  t  =  2.  The  general  behavior  shown  in  Tables  10  through 
12  is  that,  for  a  fixed  value  of  <p>  the  bound  increases  to  a  maximum  and  then 
decreases  as  the  attack  size  increases.  It  should  be  noted  that  for  <p>  as 
small  as  half  the  breakeven  value,  the  bound  never  exceeds  some  ten  percent  of 
the  population.  The  location  and  magnitude  of  the  maximum  value  of  the  bound 
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Table  10 


Upper  Bound  on  Incremental  Fatalities  -  5  psi  Hardness 


2 . 

,  2 

F(tNs5^[i4]/A*)  -  F(tNs50  J/A*)  for  ^  = 

tNs50,A/A* 

F(tNs50,A/A*> 

6 

! 

10 

30 

50 

0.05 

0.1406 

0.2358 

0.1570 

0.0431 

0.0113 

0.1 

0.2090 

0.3094 

0.2121 

i 

0.0605 

0.0160 

0.2 

0.3043 

0.3729 

0.2666 

0.0806 

0.0217 

0.3 

0.3740 

0.3937 

0.2905 

0.0920 

0.0250 

0.4 

0.4296 

0.3966 

0.3003 

0.0990 

0.0272 

0.5 

0.4760 

0.3906 

0.3024 

0.1032 

0.0287 

0.6 

0.5156 

0.3800 

0.2999 

0.1057 

0.0297 

0.7 

0.5501 

0.3669 

0.2946 

0.1069 

0.0303 

0.8 

0.5806 

0.3528 

0.2876 

0.1072 

0.0306 

0.9 

0.6077 

0.3382 

0.2795 

0.1068 

0.0307 

1.0 

0.6321 

0.3236 

0.2707 

0.1059 

0.0307 

1.5 

0.7250 

i 

0.2571 

0.2257 

0.0976 

0.0293 

2.0 

0.7872 

0.2047 

0.1857 

I 

0.0870 

0.0269 
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Upper  Bound  on  Incremental  Fatalities  -  7  psi  Hardness 


F(tNs5£[1+;]/A*)  -  F(tNs50  j/A*)for  = 


ls50,A/A* 

30 

50 

80 

0.05 

i 

0.1406 

0.2183 

0.0777 

0.0365 

0.1 

0.2884 

0.1080 

0.0514 

0.0129 

0.2 

0.3043 

0.1413 

0.0686 

0.0175 

0.3 

0.3740 

0.3730 

0.1591 

0.0786 

0.4 

0.4296 

0.3779 

0.1690 

0.0220 

0.5 

0.4760 

0.3739 

0.1743 

0.0885 

1 

0.0232 

0.6 

0.5156 

0.3652 

0.1767 

0.0240 

0.7 

0.5501 

0.3539 

0.1770 

0.0920 

0.0245 

0.8 

0.5806 

0.3413 

0.1759 

0.0924 

0.0248 

0.9  ! 

0.6077 

0.3281 

0.1739 

0.0922 

0.0249 

1.0 

0.6321 

0.3147 

1 

0.1711 

0.0916 

0.0249 

1.5 

0.7250 

0.2524 

0.1524 

0.0850 

0.0238 

2.0 

0.7872 

0.2021 

0.1320 

0.0762 

0.0219 

(as  a  function  of  attack  size)  are  given  by  Equations  292,  295  and  296  of 
Ref.  1  and  are  listed  in  Table  13  for  the  cases  considered  in  Tables  10,  11, 

12.  In  Ref.  1,  we  solve  the  inverse  problem  of  determining  the  value  of  <p> 
needed  to  guarantee  that  the  maximum  value  of  the  bound  never  exceeds  a 
specified  value;  the  results  are  shown  in  Table  14  (reproduced  from  Ref.  1) 
and  include  the  breakeven  value  ij>*  for  <p>/(F/Y).  These  results  clearly  are 
concerned  with  bounding  the  additive  effect  of  delayed  casualties.  A  comple¬ 
mentary  and  equally  informative  set  of  results  concerns  bounding  the  multi¬ 
plicative  effect  of  delayed  casualties.  We  turn  now  to  deriving  these  results. 

From  Equations  5,  6,  11  and  49,  it  follows  that 


Equation  51  clearly  is  equivalent  to  the  condition  that  the  total  combined- 
effects  casualties  from  any  mixed  attack  not  exceed  (1+n)  times  the  prompt- 

A  A 

effects  casualties  from  a  pure  airburst  attack.  For  F  =  Fq,  we  write  a  =  &0 
and  illustrate  Equation  51  for  n  =  0.5  and  Ap$  =  5,  7,  10  psi;  the  results 
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Location  and  Magnitude  of  Maximum  Value  of  Upper  Bound 
on  Incremental  Fatalities 


APS,50’  P$i 

<p> 

. F/Y 

T^S50 

Attack  Size, — , 

Producing  Maximum 

Maximum  Value 
of 

Upper  Bound 
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5.0765 
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0.0307 
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0.6694 
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50 

0.4913 

0.8154 

0.0924 

80 
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0.9474 

0.0249 

10 

10 

7.2373 

0.3106 
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30 
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50 

1.0536 
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80 

0.4738 

0.8205 

.  0.0897 

100 

0.2806 

0.8823 

0.0573 

150 

0.0229 

0.9887 

0.0053 

are  summarized  in  Table  15  and  display  the  bound  of  Equation  51  as  a  function 
of  attack  size.  It  is  apparent  that,  tor  fixed  values  of  n  and  Apg  ^q,  the 
required  protection  factor  decreases  steadily  with  increasing  attack  size. 

For  F  =  F  ,  it  is  readily  seen  that 

6  >  (l+n)1/X  -  1  (53) 


where  the  right-hand  side  of  Equation  53  is  the  limit  of  Equation  52  as  attack 
size  becomes  vanishingly  small.  It  follows  immediately  that 
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In  Table  16  we  list  this  bound  on  protection  factor  for  several  values  of 
Aps,50  and  If  the  Protection  factor  satisfies  Equation  54,  it  clearly 
satisfies  Equation  51  also;  the  converse,  of  course,  is  not  true,  as  may  be 
seen  by  comparing  Tables  15  and  16. 

From  the  results  presented  in  Tables  10  through  16,  we  see  that  the 
breakeven  values  of  protection  factor  are  modest  and,  moreover,  that  even 
1 ess-than-breakeven  values  typically  do  not  materially  alter  the  outcome, 
i.e.,  maximum  number  of  fatalities.  Consequently,  for  attacks  against  popu¬ 
lation,  the  influence  of  delayed-effects  hardness  is  generally  weak,  compared 
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Table  15 


Table  16 


to  the  influence  of  prompt-effects  hardness  and  dispersion.  For  populations 
with  protection  factor  very  small  compared  to  the  breakeven  value,  this  con¬ 
clusion  has  not  been  demonstrated  above.  If,  for  some  reason,  it  is  necessary 
to  consider  such  populations,  tighter  bounds  than  those  employed  above  may 
yet  enable  the  conclusion  to  be  extended  even  to  many  such  populations.  In 
deriving  the  bounds  used  above,  we  have  significantly  overestimated  both  the 
size  of  the  fallout  region  and  the  population  contained  therein.  Both  these 
overestimates  can  be  tightened  up,  e.g.,  by  restricting  the  number  of  effec¬ 
tively  collocated  bursts  to  be  less  than  a  small  fixed  maximum,  by  narrowing 
the  range  of  admissible  meteorological  conditions,  by  accounting  for  overlap 
between  prompt  and  delayed  effects,  by  accounting  for  extensive  unpopulated 
areas. 

The  positive  interpretation  of  our  results  concerning  delayed-effects 
hardness  is  that— precisely  because  the  breakeven  values  of  protection  factor 
are  modest— realizable,  finite  levels  of  fallout  shielding  provide  as  much 
guaranteed  protection  as  any  level  of  shielding  can  provide.  Consequently, 
efforts  to  provide  the  population  with  this  assured  protection  have  demonstrable 
payoff. 

The  bounds  derived  above  for  protection  factor  depend  upon,  inter 
alia,  "universal"  isodose  contours  and  the  area  contained  therein  (cf.  Ref.  1). 
In  Appendix  B,  we  present  numerical  details,  omitted  from  Ref.  1,  to  illustrate 
more  comprehensively  the  behavior  of  the  isodose  contours. 
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3.  CONCLUSIONS  AND  RECOMMENDATIONS 

In  this  section,  we  summarize  the  major  conclusions  which  have  emerged 
from  this  investigation  and  offer  the  corollary  recommendations  for  future 
action. 


3.1  Conclusions 

The  principal  conclusions  from  the  results  presented  here  are 
listed  below. 

1.  The  feasibility  of  interrelating  analytically  the  major  ingre¬ 
dients  of  damage/casualty  prediction  has  been  demonstrated  by  deriving  tract¬ 
able,  explicit  analytic  representations  of  the  interplay  among  the  characteristics 
of  the  attack,  targets,  environment  and  outcome. 

2.  Our  approach  for  determining  the  effective  blast  protection  of 
a  population  with  mixed,  non-uniform  protection  produces  stable,  consistent 
values  of  protection  hardness  and  accurate  estimates  of  casualties. 

3.  Relocation  schemes  which  evacuate  a  region's  population  to 
locations  where  at  least  some  people  already  reside  have  been  characterized 
abstractly  and  analyzed  parametrically,  in  particular,  the  conditions  for, 
and  consequences  of,  self-contained  regional  relocation  have  been  determined. 

4.  It  is  highly  unlikely  that  any  realistic  aporoach  could  relocate 

the  bulk  of  the  nation's  population  at  population  density  less  than  1,0C3  to 
2 

1,500  people/mi  ,  i.e.,  some  2  to  3  times  the  population  density  corresponding 
to  spreading  the  entire  population  uniformly  over  the  entire  populated  area. 
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5.  To  achieve  moderate  survivor  fractions  for  attacks  against 
population  with  several  thousand  blast-equivalent  megatons,  significant 


blast  hardening  is  required,  even  with  the  most  defense-optimistic  assumptions 
about  relocation. 

6.  For  attacks  against  population,  the  influence  of  delayed- 
effects  hardness  is  generally  weak,  compared  to  the  influence  of  prompt- 
effects  hardness  and  dispersion. 

3.2  Recommendations 

In  view  of  the  results  and  conclusions  presented  above,  we  recommend 

that 

1.  the  demographics  of  residential  population  be  analyzed  to 

a)  explore  the  validity  of  the  apparent  universality  of 
normalized,  rank-ordered  population  distribution  for  different  geographic 
regions, 

b)  provide  a  basis  for  interpreting  the  rank-ordered  popula¬ 
tion  distribution  in  other  commonly  used  formats  (e.g.,  in  terms  of  political 
subdivisions) 

2.  the  general  analytic  approach  and  results  be  extended,  refined, 
and  applied,  e.g.,  by  considering 

a)  attacks  more  general  than  attacks  against  population, 

b)  tighter  bounds  on  combined-effects  casualties, 

c)  specific  candidate  relocation  schemes. 
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APPENDIX  A 


RANK-ORDERED  POPULATION  DISTRIBUTION  FOR  SELF-CONTAINED 
RELOCATION  WITH  COMPLETE  HOSTING 

In  this  Appendix,  we  complete  the  characterization--discussed  in  Section 
2.3  of  the  text--of  the  rank-ordered  population  distribution  G  for  self- 
contained  (r=S)  relocation  with  complete  hosting  and  no  buffer  region.  The 
previous  discussion  has  shown  that  G  consists  of  four  portions,  with  pro¬ 
gressively  smaller  post-relocation  population  densities  in  each  portion.  We 
follow  the  notation  introduced  in  the  text  and  prescribe  below  the  end¬ 
points  of  each  of  these  four  portions,  i.e.,  we  state  the  incremental  area 
and  incremental  population  of  each  portion.  The  behavior  between  end-points 
is  readily  inferred  from  the  definition  of  the  relocation  and  from  the  pre¬ 
viously  discussed  properties  of  F  =  FQ  and  of  the  critical  hosting  parameters. 
The  cumulative  distribution  defined  by  the  four  portions  is  simply  G  itself. 

The  initial  portion  of  G  consists  of  that  portion  of  the  evacuation 
region  which,  even  after  relocation,  has  population  density  greater  than  u*. 
The  size  of  this  portion  is  clearly  given  by  A/A*  =  a-| »  where 

(I-y)F1 (a,)  «  F‘(a)  (1) 

Therefore  this  first  portion  of  G  has  incremental  area 
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(2) 


a(A/A*)  =  a1 


and  incremental  population 


aG  *  (l-Y)F(a1)  (3) 

The  next  portion  of  G  consists  of  that  part  of  the  host  region  which 
has  constant  population  density  u>*.  This  portion  has  incremental  area 


a(A/A*)  =  a-a 


(4) 


and  incremental  population 


aG  =  F'(a)[a-a]  (5) 

We  have  seen  that,  after  relocation,  the  smallest  population  density  is 
(1-y)w*  in  the  evacuation  region  and  is  K*F' (10)fttot/A*  in  the  host  region. 
Therefore,  we  have  from  Equation  48  of  the  text  that  the  least  densely  popu¬ 
lated  (post-relocation)  cells  are  in  the  host  region  if 

O-Y)F'(a)  >  F  ’  (10)  (6) 
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and  in  the  evacuation  region  if  the  reverse  inequality  holds.  For  a  fixed  y. 
Equation  6  is  satisfied  for  sufficiently  small  a  and  hence  for  sufficiently 
small  K*. 

If  Equation  6  is  satisfied,  then  the  portion  of  6  following  the  constant- 
population-density  portion  consists  of  the  balance  of  the  evacuation  region, 
plus  that  portion  of  the  host  region  wherein  the  post-relocation  population 
density  exceeds  (l-y)co*.  The  incremental  area  of  this  portion  is 

a 

a(A/A*)  =  a-a^+aj-a  (7) 

and  the  incremental  population  is 


AG  =  (l-y)[F(a)-F(a1)]  +  k*[F(a3)-F(£)]  (8) 

where  the  parameter  a^  is  given  (via  Equation  48  of  the  text)  by 

F *  ( a 3 )  s  (W)F'(a)  (9) 

and  thus  defines  the  host-region  location  where  the  post-relocation  population 
density  is  (l-y)u)*. 

If  Equation  6  is  satisfied,  then  the  remaining  portion  of  G  con¬ 
sists  of  the  balance  of  the  host  region  and  has  incremental  area 
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00) 


a(A/A*)  =  10-a3 

and  incremental  population 

AG  =  K*[F(10)-F(a3)]  (11) 

If  Equation  6  is  not  satisfied  (i.e.,  if  a  or  K*  is  sufficiently  large), 
then  the  portion  of  G  following  the  constant-population-density  portion  con¬ 
sists  of  the  balance  of  the  host  region  plus  that  portion  of  the  evacuation 

* 

region  wherein  the  post-relocation  population  density  exceeds  K*F' (10)ntQt/A*. 
The  incremental  area  of  this  portion  is 

a(A/A*)  =  a2-a]  +  10-a  (12) 

and  the  incremental  population  is 

AG  =  (l-Y)[F(a2)-F(a1)]  +  K*[F(10)-F(a)]  (13) 

where  the  parameter  a2  is  given  by 

(l-v)F' (a2)  =  K*F‘ (10)  (14) 

and  thus  defines  the  evacuation-region  location  where  the  post-relocation 
population  density  is  K*F' (10)ntot/A*. 
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If  Equation  6  is  not  satisfied,  then  the  remaining  portion  of  G  consists 
of  the  balance  of  the  evacuation  region  and  has  incremental  area 


a(A/A*)  =  a-a2 


and  incremental  population 


AG  =  (l-Y)[F(a)-F(a2)] 


In  sum,  then,  the  distribution  6  has  the  cumulative  behavior  determined 
by  the  four  constituent  portions.  Thus  the  first  portion  (from  Equation  2  for 
A  and  from  Equation  3  for  G)  ends  with 


A/ A*  *  a.j  G=(l-7)F(a1) 


while  the  second  portion  (from  Equations  2  and  4  for  A  and  from  Equations  3 
and  5  for  G)  ends  with 


A/A*  =  a-j+a-a  G  =  (l-'r)F(a1)  +  F'(a)[a-a] 


If  K*  is  sufficiently  small  that  Equation  6  is  satisfied,  then  the  third 
portion  (from  Equations  2,  4  and  7  for  A  and  from  Equations  3,  5  and  8  for  G) 
ends  with 
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(19) 


A/A*  =  a3  G  =  (l-y)F(a)  +  F‘ (a)[a-a]+K*[F(a3)-F(a)] 

while  the  fourth  portion  (from  Equations  2,  4,  7  and  10  for  A  and  from 
Equations  3,  5,  8  and  11  for  G)  ends  with 


A/A*  -  10  G  =  (l-y)F(a)  +  F'(a)[a-a]  +  K*[FO0)-F(a)]  (20) 


From  Equations  41  and  43  of  the  text,  we  see  that  for  the  conditions  considered 

A 

A  A 

here  (i.e.,  r=S,  A=A)  the  expression  for  G(10)  given  by  Equation  20  is  in  fact 
F(10) ,  as  it  should  be. 

A 

If  K*  is  sufficiently  large  that  Equation  6  is  not  satisfied,  then  the 
third  portion  (from  Equations  2,  4  and  12  for  A  and  from  Equations  3,  5  and 
13  for  G)  ends  with 

A/A*  =  10-(a-a2)  G  =  (l-y)F(a2)  +  F'(a)[a-a]  +  K*[F(10)-F(a)]  (21) 

while  the  fourth  portion  (from  Equations  2,  4,  12  and  15  for  A  and  from  Equations 
3,  5,  13  and  16  for  G)  ends  with 

*  *  a 

A/A*  =  10  G  s  (l-y)F(a)  +  F' (aUS-a]  +  K*[F(10)-F(a)]  (22) 

Af  before,  G(10)  as  given  by  Equation  22  is,  in  fact,  F(1 0) . 


For  y  *  1  (i.e.,  for  100  percent  evacuation)  the  distribution  6  is  de¬ 


generate  and  is  given  simply  as 


A  A 

I  [F1  («)]$*•  0  <  <  i-i 

■  -  ■  *  ,  ‘  <23) 
F'(a)[i-0  +  K*[F(a  +  fr)  -  F<a)]  J-a  <  jjr  <  10-a 

where  the  parameters  a,  a,  etc.  are  tabulated  in  Table  8  as  a  function  of  K*. 

We  now  illustrate  the  behavior  of  G  for  less  than  complete  (y=1 )  evacuation. 
Specifically,  we  consider  the  case  y  =  0.8,  corresponding  to  80  percent  evacua¬ 
tion.  For  y  =  0.8,  Equation  6  is  satisfied  for  a  less  than  5.3777  and  hence, 

A 

from  Table  9,  for  K*  less  than  a  value  intermediate  between  7  and  8.  In  Table 
A.l ,  we  list  the  values  of  a-j,  a2,  a3  determined,  for  y  -  0.8,  from  Equations 
1,  14,  9  respectively.  These  values  together  with  the  values  tabulated  for 
a,  a,  etc.  in  Table  9,  suffice  to  determine  the  distribution  G.  In  Table  A. 2 
we  summarize,  for  y  -  0.8  and  various  values  of  K*,  the  location  of,  and  the 
magnitude  of  G  at,  the  end  of  the  four  constituent  regions. 


$ 


» 


89 


Table  A. 2 


0.0321 

0.6100 

3.0368 

6.3211 


0.0850 

1.4037 

4.1597 

4.3516 


0.1337 

2.0880 

4.8193 

2.9590 


0.1747 

2.6731 

5.2710 

1.8812 


0.2088 

3.1809 

5.6069 

1.0034 


0.2372 

3.6284 

5.8698 

0.2646 


APPENDIX  B 


UNIVERSAL  ISODOSE  CONTOURS 

In  this  Appendix,  we  present  numerical  details,  omitted  from  Ref.  1, 
to  illustrate  more  comprehensively  the  behavior  of  "universal"  isodose  con¬ 
tours  as  functions  of  downwind  location,  x,  and  crosswind  location,  y. 

We  have  shown  in  Ref.  1  that,  in  the  far-field  approximation  to  WSEG-10 
fallout  phenomenology,  it  is  possible  to  characterize  each  isodose  contour 

A 

in  purely  geometric  terms  by  specifying  the  distant  location,  x,  at  which 

the  isodose  contour  crosses  the  downwind  axis.  The  crosswind  extent,  y.  .  , 

isodose 

of  the  isodose  contour  is  given  in  dimensionless  form  as 

s  yisodose  _  ,,1/2  x.  ✓  2S_  (1  -  *)  +  c  log(x)  >1/2  (1) 

"  (oHSc/W)i  x  WT  x  x 

where 

c  =  1.382  (2) 

Equation  1  presents  the  isodose  contour's  dimensionless  crosswind  extent,  v, 

A  a 

as  a  function  of  the  parameter  x/WT  and  of  the  normalized  location,  x/x.  Table  B.l 

A  A 

lists  v  as  a  function  of  x/x  for  several  values  of  x/WT.  For  some  of  these 
values  of  x/WT,  Figure  B.l  (from  Ref.  1)  presents  the  corresponding  plot  of 
v  as  a  function  of  x/x. 


Table  B.l 

Normalized  Isodose  Contours  for  Several  Values  of  x/WT  and  Normalized  Distance 


Table  B.1,  continued 


V, 

for  x/WT  = 

| 

1.0 

2.0 

4.0 

6.0 

8.0 

10.0 

0.3250 


0.4290 

0.9219 


.2377 

.5629 

.6838 

.6770 

.5745 

.3933 

.1435 

1.8306 

0.4557 


0.4725 


1 .7189 


2 


0.9044 


(x(M/°SHo)]/8SoposU 


and 


For  a  given  value  of  x/WT,  v  increases  from  0  to  a  maximum,  v  , 

then  decreases  to  0  as  x/x  increases  from  0  to  u*  (i.e.,  the  value  of  x/x 

yielding  v  )  to  1.  The  downwind  location,  x/x  =  u*,  and  magnitude,  vm  . 

md  x  iTia  x 

of  the  maximum  crosswind  extent  of  the  universal  isodose  contour  are  given  by 


x  _  (c/3)(l+2  log  u*) 


(3) 

(4) 


It  follows  immediately  from  Equations  3  and  4  that  as  x/WT  increases,  then  u* 
-1  /2 

increases  from  e  '  -  0.6065  to  2/3  =  0.6667  and  vmax  increases  monotonically 

and  unboundedly  from  (c/e)1//2  =  0.7130.  In  Table  B.2,  we  list  u*  and  vmav 

nlaX 

A 

for  several  values  of  x/WT. 

To  prepare  isodose  contours  in  the  more  conventional  format  (i.e.,  as 
loci  of  points  in  the  x-y  plane),  it  suffices  to  recognize  that 


x  _  x  x 
WT  '  WT  “ 


(5) 


^isodose 

WT 


(oHSc/W>  WT  v 


(6) 


Table  B.3  lists  values  of  y^S0(j0Se/wT  as  a  function  of  x/WT  for  o^Sc/K  =  0.1 
and  for  several  values  of  x/WT.  Figures  B.2  and  B.3  (from  Ref.  1)  present  the 
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Table  B.3 


Universal  Isodose  Contours  for  o^Sc/W=0.1  and  Several  Values  of  x/WT 


x/WT= 

-0.1 

x/WT=0.2 

x/WT=0.4 

x/WT=0.6 

x/WT 

y/WT 

x/WT 

y/WT 

x/WT 

y/WT 

0.099 

0.0017 

0.198 

0.0035 

0.396 

0.0075 

0.594 

0.0118 

0.095 

0.0037 

0.19 

0.0076 

0.38 

0.0162 

0.57 

0.0256 

0.09 

0.0050 

0.18 

0.0104 

0.36 

0.0219 

0.54 

0.0346 

0.08 

0.0065 

0.16 

0.0134 

0.32 

0.0282 

0.48 

0.0444 

0.07 

0.0072 

0.14 

0.0147 

0.28 

0.0310 

0.42 

0. 0487 

0.06 

0.0073 

0.12 

0.0150 

0.24 

0.0316 

0.36 

0.0495 

0.05 

0.0071 

0.10 

0.0145 

0.20 

0.0304 

0.30 

0.0476 

0.04 

0.0065 

0.08 

0.0133 

0.16 

0.0278 

0.24 

0.0433 

0.03 

0.0056 

0.06 

0.0114 

0.12 

0.0237 

0.18 

0.0367 

0.02 

0.0043 

0.04 

0.0087 

0.08 

0.0180 

0.12 

0.0279 

0.01 

0.0026 

0.02 

0.0052 

0.04 

0.0106 

0.06 

0.0164 

0.005 

0.0015 

0.01 

0.0029 

0.02 

0.0060 

0.03 

0.0092 

0.001 

0.0004 

0.002 

0.0007 

0.004 

0.0015 

0.006 

0.0022 

Figure  B.3  Isodose  Contours  for  a  HSr/W  *  0.1  and  x/WT  >  1 


corresponding  plots.  To  obtain  results  for  any  other  value  of  ouS  /W,  it 

H  C 

suffices  to  multiply  the  scale  of  the  crosswind  axis  by  (o  S  /W)/0.1.  It 

H  C 

is  clear  that  the  downwind  location  and  magnitude  of  the  maximum  crosswind 
extent  of  each  isodose  contour  are  given  by 


JL  * 
WT 


WT  "  WT 


<°hVw>  St  *, 


max 


(7) 

(8) 
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